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ABSTRACT 


An  experimental  study  of  the  development  of  air 
velocity  profiles  in  a  vertical  6.2  55  inside  diameter  aluminum 
tube  with  a  diverging  conical  entrance  was  undertaken.  Fully 
developed  velocity  profiles  were  obtained  with  a  total  head 
pitot  tube  for  Reynolds  numbers  ranging  from  9000  to  56000. 

The  development  of  the  center-line  velocity  along  the  tube 
axis  was  measured  with  a  travelling  carriage. 

Various  means  for  correlating  velocity  profiles 
were  investigated,  from  the  simple  power  lav/,  universal 
velocity  distribution,  Gill  and  Scher  '  s  improved  correlation 
(based  on  a  modified  mixing  length  theory)  to  the  correlation 
derived  by  Pai .  The  different  relations  which  have  been  de¬ 
rived  for  the  region  near  the  wall  were  analyzed  and  compared. 

The  experimental  results  correlated  very  well  with 
an  equation  of  the  form  suggested  by  Bogue,  a  modified  uni¬ 
versal  velocity  distribution  which  takes  into  account  devia¬ 
tions  from  the  universal  velocity  distribution  with  radial 
position.  However,  no  significant  effect  of  Reynolds  number 
could  be  observed  over  the  limited  range  investigated.  The 
equation  derived  by  Pai  also  correlated  the  data  very  well  and 
since  it  has  the  advantage  of  being  valid  from  the  center- line 
to  the  wall,  it  would  seem  to  be  the  most  advantageous  cor¬ 
relation  over  the  range  of  Reynolds  numbers  investigated. 


The  center-line  velocity  was  found  to  reach  its 
maximum,  and  constant  value,  at  40  diameters  downstream  from 
the  entrance.  The  hydrodynamic  starting  length  v/as  found  to 
increase  with  increasing  Reynolds  number. 
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I. 


INTRODUCTION 


Experimental  studies  on  the  development  of  velocity 
profiles  for  air  flowing  upwards  in  a  vertical  tube  with  a 
diverging  conical  entrance  were  carried  out  by  Bowers(l)  and 
Sandall(2)  in  the  Department  of  Chemical  and  Petroleum  Engineer¬ 
ing  at  the  University  of  Alberta.  Their  studies  encompassed 
air  flows  ranging  in  Reynolds  numbers  from  roughly  1000  to 
30,000.  However,  both  of  these  earlier  workers  were,  to  some 
extent,  limited  in  the  scope  and  generality  of  their  measure¬ 
ments  . 

The  bulk  of  the  difficulties  which  they  encountered 
originated  with  an  air  supply  subject  to  excessive  ambient 
temperature  and  pressure  changes,  and  from  the  lack  of  a 
suitable  means  of  calibrating  their  hot  wire  anemometer  at  low 
velocities.  To  minimize  these  difficulties,  the  following 
changes  have  been  made  to  the  experimental  equipment: 

1.  an  air  compressor  of  higher  compression  and 
volumetric  capacity,  Re  up  to  55,000  in  the  test  section, 
has  been  installed; 

2 .  with  an  after-cooler  the  air  supply  temperature 
could  be  controlled  within  ~1°F; 

3.  further  changes  in  the  vertical  tube  entrance 
geometry  were  necessary  to  assure  radial  symmetry  of  the 
developing  velocity  profiles; 
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4.  calibration  of  secondary  velocity-measuring 
devices  was  eliminated  by  using  a  pitot  total  head  tube  of 
very  small  diameter  for  all  air  velocity  measurements  (Re  from 
9,000  to  55, 000) ; 

5.  a  travelling  probe  which  measured  center-line 
velocity  was  used  to  indicate  the  velocity  profile  development 
along  the  tube  axis. 

Using  the  improved  experimental  equipment,  the  fol¬ 
lowing  objectives  were  undertaken  in  this  study: 

1.  to  determine  the  velocity  profiles  for  various 
values  of  Reynolds  numbers  corresponding  to  experimentally 
verified  fully  developed  turbulent  flows; 

2 .  to  define  the  distances  along  the  tube  axis  in 
which  fully  developed  flows  would  be  encountered; 

3.  to  review  and  test  various  published  velocity 
correlations  for  their  suitability  in  correlating  the  measured 
fully  developed  profiles; 

4.  to  examine  the  axial  development  of  the  velocity 
profile  in  the  entrance  region  using  constancy  of  the  center- 
line  velocity  as  the  criterion  rather  than  constancy  of  the 
wall  shear  stress. 

The  problem  of  correlating  experimentally  measured 
velocity  profiles  in  generalized  mathematical  terms  has  been 
and  continues  to  remain  the  subject  of  considerable  discussion. 
A  truly  universal  function  which  relates  local  mean  velocity 
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to  radial  position  is  not  available  hence  the  worker  in  this 
area,  faced  with  the  need  for  a  velocity  profile  relationship, 
is  restricted  to  using,  in  direct  or  modified  form,  the 
various  relationships  which  have  been  proposed.  The  following 
section  discusses  ways  of  correlating  velocity  profiles  in 
tubes  with  special  reference  to  the  flow  of  Newtonian  fluids 
such  as  air . 
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II.  THEORY 

1.  The  Correlation  of  Experimental 

Velocity  Profiles  in  Circular  Ducts 

In  the  case  of  laminar  flow,  theoretical  relation¬ 
ships  to  predict  velocity  profiles  have  been  derived.  Tur¬ 
bulent  flow  cannot  be  analyzed  satisfactorily  from  a 
theoretical  standpoint  and  so,  only  empirical  or  semi-empirical 
relationships  to  predict  velocity  profiles  in  fully  developed 
turbulent  flow  have  been  obtained.  The  condition  of  fully 
developed  flow  is  referred  to  as  the  condition  where  the 
velocity  profile,  and  in  particular  the  center-line  velocity, 
has  reached  its  "constant"  value.  The  various  approaches 
to  correlate  mean  point  velocity  in  turbulent  flow  are  dis¬ 
cussed  in  four  sections  as  follows: 

a.  Power  Law  Distributions 

b.  Logarithmic  Distributions  Including  the  Universal 
Velocity  Distribution 

c.  Improved  Universal  Velocity  Distributions  Near 
the  Wall 

d.  Other  Distributions  in  the  Turbulent  Core 
a .  Power  Law  Distributions 

The  power  law  of  velocity  distribution  for  turbulent 
flow  in  circular  ducts  represents  a  simple  way  of  analyzing 
data.  Using  a  power  function  with  a  single  term,  the  power 


law  is  given  by 


u 

u 


m 


(~)n 

R 


(1) 


This  form  was  tested  by  Nikuradse(6)  and  was  found 
to  correlate  his  experimental  data  well,  however,  the  exponent 
n  varied  from  6  at  a  Reynolds  number  of  4000  to  10  at  a  Rey¬ 
nolds  number  of  3,240,000. 

It  is  interesting  to  note  that  equation  (1)  can  be 
derived  from  the  Blasius  friction-factor  equation  which  relates 
the  friction  factor  to  the  Reynolds  number  in  the  following 
way 


f  =  C  (Re) 


■1/4 


(2) 


Using  this  approach,  Prandtl(3)  obtained 

1 

u  y  7 

—  =  C'  (-)  (3) 

u  R 

m 

which  is  only  valid  for  Reynolds  number  up  to  100,000  for 

which  the  Blasius  friction-factor  equation  is  applicable. 

With  a  reliable  functional  relationship  between  the 

exponent  n  and  the  Reynolds  number  the  power  law  would 

represent  a  useful  and  simple  way  of  correlating  velocity  data. 

b .  Logarithmic  Distributions  Including  the  Universal 
Velocity  Distribution 

Several  velocity  distribution  equations  have  been 
derived(3)  in  the  form  of  semi-logarithmic  functions.  One 


■ 
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of  the  earliest  approaches  was  that  of  Prandtl.  Using  his 
mixing  length  theory,  he  defined  the  total  shear  stress  by 
the  sum  of  the  viscous  and  the  turbulent  shear  stresses  as 
follows: 


(1  du  p  du  ~ 

t  =  t  +t  =  —  —  +  —  (1  — )  (4) 

Uc  dy  gc  dy 

Neglecting  viscous  shear,  and  assuming  the  shear  stress  to  be 
constant  and  equal  to  the  shear  stress  at  the  wall,  Prandtl 
obtained 


P  du  2 

,  =  XW  =  -  (1-) 

dy 


By  introducing  the  friction  velocity,  u* 
equation  (5)  can  be  simplified  to 

du 


(5) 


=  ( 


lw 


!i)1/2 


u*  =  1 


dy 


(6) 


If  the  mixing  length,  1,  is  assumed  to  be  a  linear  function 
of  y,  equation  (6)  may  be  integrated,  yielding 

u* 

u  =  —  In  y  +  C  (7) 

K 

Prandtl  employed  the  boundary  condition,  u  =  u  ,  for  y  =  R 
and  obtained 


y 

u  =  u_  +  k  u*  In  —  (8) 

R 

In  spite  of  the  number  of  drastic  assumptions  which  were 
incorporated, its 
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predictions ( 3)  compare  well  with  the  experimental  data. 


von  Karman(3)  showed  that  for  flow  at  a  distance 


from  the  wall, 


the  shear  stress  is  given  by 
2  du  A 


T  = 


d2u  2 

<ry2> 


(9) 


Introducing  Prandtl's  mixing  length  theory,  and  rearranging, 
gives 


I< 


du 

dy 


d2U 


dy 


(10) 


y, 


By  using  the  linear  shear  stress  distribution,  t  =  (1  -  —)  tw, 

R  w 


in  equation  (9)  and  rearranging,  the  following  differential 
equation  was  obtained. 


—  -  \f ^ 

dy2  V  gc  tw 

After  integrating  equation  (11)  twice,  and  applying  the  boundary 
conditions, 

oo  ,  for  y  -  0 

u  ,  for  y  =  R, 


du 

dy 

and 

u  = 


1  -i 


du  9 

k  (— r 

dy 


(11) 
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von  Karman  obtained, 


1 

U  =  U„  +  —  u* 

In  (1  -  * 

\jl  -  -  ) 

+  M 

Ci' 

m  K 

V  R 

V 

R 

The  constant,  K,  appearing  in  the  equations  of  both  Prandtl 
and  von  Karman,  is  considered  to  be  a  universal  constant  and 
its  value  has  been  experimentally  determined  at  small  distances 
from  the  wall (6) . 

Wang(3),  using  a  similar  approach  involving  an  ex¬ 
pression  for  the  mixing  length  based  on  an  actual  velocity 
distribution,  obtained 


u  = 


where 


f(y,R)  = 


(13) 


2  tan 


1  - 


y 

R 


y  a  y 

1  -  -  +  1 . 7  5  \/  1  -  -  +  1.53 

R  V  R 


-  0.572  In 


y  A  |  y 

1 - 1.75\l  -  -  +  1,53 

r  Nr 


+ 


1.14  tan 


1.53  -  (1  -  -) 

R 


(14) 
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Equation  (13)  supposedly  agrees  better  with  experimental  data, 
however,  the  second  term,  f(y,R),  is  very  complex  and  incon¬ 
venient  to  use . 

A  limited  comparison  of  the  three  logarithmic 
velocity  distribution  equations  (8),  (12)  and  (13)  is  presented 

by  Knudsen  and  Katz (3) .  The  few  data  which  were  plotted  in¬ 
dicate  good  agreement  with  equation  (8)  at  low  Reynolds  number 
and  very  good  agreement  with  equation  (13)  at  high  Reynolds 
number .  It  is  important  to  note  that  all  three  equations  become 
invalid  when  y  approaches  zero. 

A  more  general  correlation,  the  so-called  universal 
velocity  distribution,  employs  the  following  dimensionless 
quantities , 


+  .  + 

Differentiating  u  with  respect  to  y  yields: 

du  du'r  9c  Tw  du'1' 

- —  =  u*  — —  =  - —  • — ~ 

-J- 

dy  dy  p  dy 


2  + 


u 


du 


c  *w 


C  “W 


,2  + 
d  u 


dy 


+‘ 


(15) 


(16) 


(17) 


(18) 


"S?  dy'dy 


2 

P 


P 


» 


10  - 


Combining  equations  (17)  and  (18)  with  equation  (9),  in  which 
viscous  shear  has  been  neglected,  and  integrating  the  resulting 
second  order  differential  equation  yields 


dy 


+ 


du 


+ 


Ky 


+ 


(19) 


Near  the  wall,  the  velocity  gradient  is  very  large  and  hence 


dy 


+ 


+ 


du 


+ 


becomes  very  small,  giving  C  =  0  for  y'  =0.  Integrating 


a  second  time  gives 


u 


K 


In  y 


(20) 


which  is  called  the  universal  velocity  distribution  equation 
for  turbulent  flow.  The  constant  of  integration,  C',  must  be 
determined  experimentally. 

Equation  (20)  is  only  valid  at  a  distance  from  the 
wall,  hence,  near  the  wall  where  the  variation  in  shear  stress 
can  be  neglected,  a  laminar  sublayer  is  assumed  to  exist.  In 
this  region,  where  viscous  shear  is  important,  Newton's  law  of 
viscosity  predicts 


-g 


T 

C  W 


du  o  du+ 

h  — ■  =  pu*  — T 

dy  dy4 


du 


dy 


+ 


and  thus,  u+  =  yf  +  C  .  Since  u+  =  0  for  yh  =  0,  the  inte¬ 


gration  constant,  C,  vanishes, giving 
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Using  the  experimental  results  of  Nikuradse ( 3) ,  equation  (20) 

was  found  to  be  valid  for  y+  p*  30.  The  limited  data  available 

for  the  regions  near  to  the  wall  suggest  that  equation  (21) 

,  + 

might  be  applicable  for  y  <  5.  In  the  region  between  this 
laminar  sublayer  and  the  turbulent  core,  the  velocity  distribu¬ 
tion  has  been  represented  by  an  empirical  equation  of  the  form 
of  equation  (20). 

The  universal  velocity  distribution  in  a  pipe  can 
be  represented  by  the  following  three  equations ( 3) , 


u+ 

=  y+ 

+ 

A 

d 

(21) 

+ 

u 

=  5.0 

In 

y+  -  3.05 

,  5  <  y+<30 

(22) 

u+ 

=  2.5 

In 

y+  +  5.5 

30  <  y+ 

(23) 

The  constants  as  well  as  the  applicable  range  were  obtained 
using  the  experimental  results  of  Nikuradse. 

This  set  of  equations,  (21)  to  (23),  is  often  used 
to  predict  the  complete  velocity  profiles  for  fully  developed 
turbulent  flows  in  pipes  or  tubes. 

In  a  more  recent  work  by  Deissler( 7 ), slightly  dif¬ 
ferent  constants  were  obtained.  For  air  flows  corresponding 
to  Reynolds  number  of  8,000  to  220,000,  Deissler  obtained  in 
the  turbulent  core,  y  >26, 

u4  =  2,78  In  y'  +  3.8 


(24) 
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However,  it  should  be  noted  that  Nikuradse's  data  cover  a  much 
wider  range  of  Reynolds  numbers,  from  4,000  to  3,240,000. 

Although  the  set  of  universal  velocity  distributions, 
equations  (21),  (22),  (23)  may  be  inconsistent  with  the 

physical  situation,  this  approach  yields  a  surprisingly  good 
correlation  of  experimental  data. 

There  is  strong  evidence (6)  that  the  flow  becomes 
laminar  very  close  to  the  wall  and  that  a  separate  turbulent 
central  core  exists.  However,  a  gradual  change  from  the  one 
to  the  other  would  be  anticipated.  Therefore,  the  concept  of 
three  distinct  zones  would  appear  to  be  incorrect.  The  derived 
equation  (20)  also  does  not  predict  a  zero  velocity  gradient 
at  the  center  of  the  pipe.  Finally,  consistent  deviations  of 
the  experimental  results  from  equation  (20)  suggests  that  u+ 
and  y+  cannot  be  the  only  dimensionless  groups  belonging  in 
the  universal  relationship. 

c .  Improved  Universal  Velocity  Distributions  Near  the  Wall 

In  the  region  near  to  the  wall,  other  correlations 
which  appear  to  be  equally  satisfactory  have  been  developed. 
Deissler(7)  derived  an  equation  for  flow  in  the  vicinity  of 
the  wall  in  which  the  effects  of  the  magnitudes  of  u  and  y 
upon  the  turbulent  shear  have  been  taken  into  account. 

The  total  shear  stress  defined  as  the  sum  of  the 
viscous  and  the  turbulent  shear  stresses  is  given  by 


+ 


P 

9c 


T 


\i  du 

gc  dy 


du 

6  - 

dy 


(25) 
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As  an  approximation,  the  following  functional  relation  for  £ 
the  coefficient  of  eddy  diffusivity,  was  assumed: 


€  =  f  (u,  y) 

From  dimensional  analysis  Deissler  obtained 

2 

t  =  n  uy 

where  n  is  a  constant  of  proportionality. 
Equation  (25)  becomes 

*  p  P  2  du 

t  =  (  +  —  n  uy)  ~ 

gc  9C  <Jy 


(26) 


(27) 


(28) 


The  dimensionless  quantities,  as  defined  by  equations  (15)  to 
(17),  substituted  in  equation  (28),  with  shear  stress  approxi 
mated  over  the  entire  sublayer  by  the  wall  shear  stress,  t  , 
gives  the  following  identity 


1  = 


/  n  2  +  +» 

( 1  +  n  u  y  ) 


du 


dy‘ 


(29) 


Upon  separating  variables  and  integrating  from  the  lower 
boundary  condition,  u+  =  0,  when  y+  =  0,  Deissler  obtained 
the  solution, 

(nu+) 2  u+  ( nu+ ) 2 

du+  (30) 

o 


Equation  (30)  can  be  expressed  in  terms  of  the  normal  error 
function  of  nu'1'  and  the  integral  of  the  normal  error  function 
(both  available  in  mathematical  tables)  to  give 
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nu+ 


(nu+) 2 

d ( nu ^ ) 

(nu+) ^ 


(31) 


1  2 
— -  e 

V 2?r 

Equation  (31)  has  been  found  to  fit  the  experimental  data  in 

.t. 

the  region,  y  26. 

Wasan  et  al(12)  derived  an  equation  for  flow  in  the 
vicinity  of  a  pipe  wall  based  upon  the  equations  of  continuity 
and  motion.  The  mean  and  fluctuating  velocity  components  were 
expanded  separately  in  Taylor  series  as  functions  of  radial 
position  and  the  equation  of  motion  was  integrated.  It  was 
shown  that  the  velocity  distribution  function,  in  dimensionless 
form,  can  be  expressed  as 


u 


=  Y 


+ 


-U  4.  A  4.  -1-  S 

U/,  (y  )  +  U,  (y  ) 


(32) 


When  the  shear  stress  is  assumed  constant  as  Deissler  did  in 

-f-  ~f * 

deriving  equation  (31),  it  was  found  that  U^,  and  are  uni¬ 
versal  constants.  To  determine  the  constants  and  the  value 
of  y  at  which  a  smooth  and  continuous  transition  to  the 
logarithmic  distribution  occurs,  the  values  of  u+  and  the 
first  and  second  derivative  of  u  v/ith  respect  to  y  ,  obtained 
from  equation  (32),  are  matched  with  corresponding  values 
obtained  from  equation  (23) .  The  transition  was  found  to 
occur  at  y*!*  =  20  and  the  velocity  distribution  for  y’r-<20  is  as 
follows 


u 


+ 


y 


+ 


1.04  (10  4)(y'!‘)4  +  3.03(10  6)(y*)b  (33) 
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Both  equation  (31)  and  (33)  approach  u  =  y  for 
y +  <  5  and  so,  they  satisfy  the  condition  of  laminar  flow 
close  to  the  wall.  Equation  (33)  was  obtained  matching 

equation  (32)  and  (2  3)  at  y'1'  =20.  In  this  region  equation  (2  3) 

*  •  .  + 

is  not  valid  and  it  can  be  expected  that  u  obtained  from 

equation  (33)  will  be  too  high. 

d .  Other  Distributions  in  the  Turbulent  Core 

To  obtain  equation  (22),  a  constant  shear  stress  was 

assumed  across  the  pipe.  Using  the  approach  of  von  Karman,  an 

alternative  equation  of  the  form  of  equation  (12)  was  derived 

by  Deissler(7) . 


u+  - 

1 

K 

In  ( 1  - 

V 

_  xl) 

R+ 

+  \ 

1 

1 

y+ 

R+ 

+ 

um 

u* 

Since  equation  ( 34)  is  expressed  in  terms  of  the  maximum 
velocity,  u  ,  it  is  no  longer  a  universal  correlation.  Further¬ 
more,  it  can  be  seen  that  equations  (12)  and  (34)  essentially 
are  identical . 

Rothfus  et  al(8,9)  from  an  examination  of  the 
available  published  data  found  that  the  relationship  between 
ur  and  y+  is  dependent  on  the  Reynolds  number.  Using  this 
concept,  they  proposed  a  modification  of  the  relationship 
between  u+  and  y+  which  would  result  in  a  unique  correlation 
of  the  steady  component  of  the  velocity.  By  examining  u' ,  y+ 
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data,  the  necessary  correction  was  made  by  multiplying  one  or 
both  of  the  coordinates  by  a  factor  which  would  be  unity  at 
high  Reynolds  numbers  and  which  would  progressively  depart  from 
unity  when  the  Reynolds  number  is  decreased.  This  correction, 
the  ratio  of  the  bulk  velocity  to  the  maximum  velocity,  was 
applied  and  it  was  found  that  the  modified  coordinates, 


u 


av 


u 


m 


and 


+ 

Y 


u 


m 


u 


av 


are  represented  by  a  single  line  over 


the  entire  turbulent  range  tested,  3000  <  Re  <  3,240,000. 
However,  in  the  vicinity  of  the  wall  further  experimental  work 
was  reported  to  be  necessary.  Discrepancies  can  be  expected 
in  the  laminar  sublayer  in  particular  if  a  true  laminar  sub¬ 
layer  exists. 

The  semi-logarithmic  universal  correlation  given  by 
equation  (20)  was  extensively  tested  by  Bogue  and  Metzner(lO). 
They  found  that  the  data  showed  small  but  consistent  deviations 
from  the  correlation.  The  deviations  are  primarily  a  function 
of  radial  position  and  to  a  lesser  extent,  of  the  friction 
factor  (or  more  generally  speaking,  of  the  Reynolds  number) . 

A  correlation  of  the  following  form  was  used, 


u '  -  C(y, f )  = 

The  correction  term,  C(y,f) 
is  given  by  the  equation 

C ( y , f )  =  0.05 


2.42  In  y+  +5.57  (35) 

was  determined  by  experiment  and 

Y  ^ 

-(-  -  0.8)2 

R 

exp  ( — - — - — )  (35) 

0.15 
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in  which  f  =  Fanning  friction  factor.  Similar  correlations 
have  been  suggested  by  Millikan ( 11) ,  Reichardt ( 11)  and  Hinze(5), 
however,  their  correction  terms  only  took  into  account  the 
deviation  as  a  function  of  radial  position,  and  assumed  no 
functional  relationship  with  Reynolds  number.  The  different 
correction  functions  are  compared  graphically  by  Bogue(ll). 
Based  on  the  data  of  Laufer  and  Nikuradse,  significant  improve¬ 
ments  in  the  correlation  using  equation  (35)  were  reported. 

Gill  and  Scher(16)  derived  a  velocity  distribution 


based  upon  a  modification  of  the  Prandtl  mixing  length  expres¬ 
sion.  Instead  of  assuming  that  1  =  ky,  the  following  mixing 
length  expression  was  assumed, 


1  = 

Ky  (1 

,  y 

-  e  ) 

(37) 

where , 

0  = 

+ 

Vm  " 

60 

(38) 
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Expression  (37)  purports  to  take  into  account  the  wall 
dampening  effect  upon  the  mixing  length.  For  the  case,  y 
approaching  R,  equation  (37)  reduces  to  1  =  Ky . 

Starting  from  equation  (4),  applying  a  linear 
variation  of  shear  stress  with  position,  and  then  using 
equation  (37)  for  the  mixing  length,  Gill  and  Scher  showed 
that 

+  .  r - 

1  +  4cd 

— -  dy 


o 


2c 


(39) 


18 


where , 


and 
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kV 


(1 
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.+ 
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d  =  1 


(40) 


(41) 


Equation  (39)  describes  the  velocity  distribution  from  the  wall 
to  the  center  of  the  duct.  A  zero  velocity  gradient  is  obtained 
at  the  center  of  the  tube  and  the  Reynolds  number  effect  upon 
the  velocity  distribution  has  been  taken  into  account.  Equation 
(39)  still  does  not  represent  a  true  universal  relationship  in 
terms  of  the  dimensionless  groups,  u~h  and  y+,  and  graphically, 
the  velocity  distribution  must  be  represented  by  a  different 
curve  for  each  value  of  the  Reynolds  number.  Equation  (39) 
has  not  been  tested  extensively  for  pipeflow,  but  it  offers 
the  great  advantage  of  being  valid  from  the  wall  to  the  center 
line  . 

Pai(13)  developed  a  correlation  for  turbulent  flow 
between  parallel  plates  that  is  valid  over  the  entire  distance 
between  the  plates.  It  has  subsequently  been  shown  by  Brodkey 
et  al(14,15)  that  the  correlation  could  be  extended  to  flows 
in  smooth  tubes,  to  Newtonian  as  well  as  non-Newtonian  fluids, 
and  it  would  still  be  valid  over  the  entire  cross-section. 

For  the  flow  of  Newtonian  fluids  in  a  tube,  the  velocity  dis¬ 
tribution  is  represented  by  a  three-term  power  series  of  the 
form 


r  2 

al  (R> 


+ 


a2 


2M 


(42) 
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Integration  of  equation  (42)  yields  an  expression  for  the  ratio 
of  average  to  maximum  velocity 


uav  1  1 

—  1  +  —  a  -|  +  — — " — —  a  p 

u  2  1  l  +  m  1 


(43) 


In  equation  (42)  and  (43) 


S  -  M 


M  -  1 


(44a) 


1 


M  -  1 


( 44b) 


By  definition 
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2u 
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(45) 


With  the  aid  of  the  following  relation 


y 


m 


Re 


64 


(46) 


S  was  computed  in  terms  of  Reynolds  number, 


S  =  0.585  +  3.172(10  3)  Re°'833 


(47) 


Equation  (47)  was  found  to  be  valid  for  Re  :>  2,800.  The 
variable  M  was  obtained  from  equation  (43),  by  trial  and 
error  solution,  and  correlated  in  terms  of  Reynolds  number 


M  =  -0.617  +  8.211  (10“3)  Re0,/86  (48) 


Equation  (48)  was  found  to  be  valid  for  all  Reynolds  numbers. 
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The  agreement  with  experimental  data  was  evidently 
good  but  certain  discrepancies  at  high  Reynolds  numbers, 

Re  100,000  were  reported. 

2  .  The  Development  of  Velocity  Profiles  in  the 
Entrance  Section  o f _ a  Cir cular  Duct 

For  laminar  flow,  relationships  to  predict  entrance 
lengths  and  velocity  profiles  in  the  entrance  section  have  been 
developed( 3) .  These  relationships  agree  well  with  experimental 
data.  For  turbulent  flow,  however,  no  reliable  relationship 
is  known.  The  entrance  length  could  basically  be  defined 
in  two  ways:  Most  commonly, it  is  defined  as  the  distance  re¬ 
quired  for  the  complete  velocity  profile  to  reach  its  final 
value;  or  it  could  be  taken  as  the  distance  required  for  the 
pressure  gradient  along  the  pipe  to  become  uniform.  In  the 
present  report,  the  first  definition  is  applied  unless  other¬ 
wise  stated. 

Latzko(18)  derived  a  relationship,  based  on  a 
rounded  entrance  with  an  initially  flat  profile,  to  predict 
the  entrance  length  as  well  as  velocity  profiles  in  the 
developing  section  for  turbulent  flow.  However,  the  results 
were  not  satisfactory,  with  the  predicted  entrance  lengths 
being  5  times  shorter  than  what  had  been  experimentally  ob- 
tained(18).  Using  a  rounded  entrance,  Nikuradse(4)  obtained, 
for  Reynolds  numbers  ranging  from  4000  to  3,000,000,  fully 
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developed  flow  at  25  to  40  diameters  from  the  entrance. 
Deissler(7),  also  working  with  a  rounded  entrance,  obtained 
fully  developed  flow  at  a  location  between  46  and  100  diameters 
for  a  Reynolds  number  of  46,000.  Kirsten(18),  working  in  the 
range  of  Reynolds  numbers  from  20,000  to  100,000,  reported  an 
entrance  length  of  40  diameters  for  a  smooth  entrance  and  50 
diameters  for  a  sharp-edged  entrance.  Furthermore,  he  found 
that  the  entrance  length  decreased  for  a  smooth  entrance  and 
increased  for  a  sharp-edged  entrance  with  increasing  Reynolds 
number.  However,  Kirsten  correlated  uav/u  versus  Reynolds 
number  and  found  considerable  scatter  in  his  results,  presum¬ 
ably  from  the  error  associated  with  obtaining  the  average 
velocity.  Kirsten  also  compared  velocity  profiles  in  the 
entrance  section  with  the  Power  Law  Distribution.  For  a 
rounded  entrance,  the  initial  profiles  were  laminar,  with  an 
exponent  close  to  1,  while  with  a  sharp-edged  entrance,  the 
initial  profiles  had  an  exponent  higher  than  that  of  fully 
developed  flo w.  It  should  be  noted  that  Kirsten  obtained 
his  flow  by  using  an  exhaust  fan  at  the  discharge  end  and 
thus  had  an  entering  stream  of  low  turbulence  level. 

Deissler  found,  in  contradiction  to  Kirsten,  that 
the  development  of  the  velocity  distribution  was  more  rapid 
for  the  sharp-edged  entrance  than  for  the  rounded  entrance. 

With  the  sharp-edged  entrance,  fully  developed  flow  was  ob¬ 
tained  after  roughly  45  diameters  and  with  the  rounded  entrance, 
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it  was  reported  that  the  distribution  was  still  developing 
at  100  diameters  from  the  entrance.  The  difference  between 
the  entrance  lengths  experienced  by  Kirsten  and  Deissler,  both 
working  with  rounded  entrances,  indicate  a  strong  effect  of 
the  turbulence  level  in  the  entering  stream  upon  the  entrance 
length.  Based  upon  the  scarce  experimental  results  mentioned, 
it  seems  apparent  that  there  are  at  least  three  conditions 
that  will  determine  the  entrance  length  for  turbulent  flow  in 
a  smooth  pipe? 

1.  entrance  configuration 

2 .  previous  history  of  the  entering  stream  or  the 
turbulence  level 

3.  Reynolds  number  based  on  the  cross-section  of 

the  pipe . 

The  distance  required  for  the  local  pressure 
gradient  along  the  pipe  to  become  steady  is  considerably  less 
than  that  required  for  the  development  of  the  complete  velocity 
profile.  Pascucci(17)  found  that  the  development  of  the 
friction  factor  or  shear  stress  is  completed  at 

-  =  3.80  log10  Re  -  2.14  (49) 

D 

Latzko  predicted! 3) 

L  1/4 

=  0.623  Re  '  (50) 

D 
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For  a  Reynolds  number  of  50,000,  equations  (49)  and  (50)  give 
an  entrance  length  of  15.7  and  9.3  diameters,  respectively, 
indicating  that  the  pressure  gradient,  and  thus  the  velocity 
profile  adjacent  to  the  wall,  becomes  established  long  before 
the  velocity  profile  in  the  core.  The  latter  was  mentioned 
earlier  in  this  report  to  be  established  at  40  to  100  diameters. 


~\ 
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III.  EXPERIMENTAL  PROGRAM 

The  equipment  used  in  this  work  has  been  described( 1 , 2 )  . 
However,  a  considerable  number  of  improvements  and  additions 
were  incorporated.  These  modifications  were  listed  briefly  in 
the  Introduction  section  and  the  present  physical  setup  is 
shown  schematically  in  Figure  1.  A  further  description  of  the 
equipment  and  its  operation  follows  in  the  Equipment  portion 
of  this  section. 

Basically,  this  equipment  provided  a  means  of 
regulating  the  flow  rate  and  the  temperature  of  air  with 
mean  flows  ranging  from  Reynolds  numbers  of  9,000  to  55,000 
(based  upon  flow  in  the  test  section) .  Fixed  locations  along 
the  pipe  axis  are  available  for  measuring  velocity  profiles. 

The  development  of  the  velocity  profiles,  a  matter  of  some 
importance  in  this  work,  was  studied  by  means  of  a  moving 
carriage  placed  within  the  tube.  The  experimental  test  pro¬ 
gram  is  described  immediately  following  the  equipment  section. 

1 .  Equipment  and  Its  Operation 
a  .  The  Test  Section 

Air  flows  were  studied  within  a  vertically  mounted 
aluminum  tube  28.5  feet  long,  6.255  inches  in  average  inside 
diameter,  and  0.203  inches  in  wall  thickness.  The  tube  was 
tapped  at  five  fixed  locations  for  velocity  profile  measure¬ 
ments  as  shown  in  Figure  1.  The  entrance  and  exit  cones 
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diverged  and  converged,  respectively,  from  a  one  inch  diameter 
entrance  to  the  6.255  inch  tube  diameter  and  then  to  a  two  inch 
diameter  exit. 

Although  the  profiles  measured  at  41.5  pipe  diameters 
down  stream  from  the  entrance  were  found  to  be  symmetrical, 
considerable  difficulty  was  experienced  in  designing  an  entrance 
section  which  produced  symmetrical  profiles  at  the  tube  inlet. 

To  eliminate  the  distortion  of  the  profile  which  results  from 
the  right  angle  tee  immediately  preceding  the  inlet  cone,  corner 
vanes  were  designed  and  installed( 2 1 ,  22) .  These  vanes  partially 
remedied  the  flow  distortion  and  reduced  the  pressure  drop 
associated  with  the  turn  but,  as  can  be  seen  in  Figure  2,  they 
were  not  completely  successful.  The  final  arrangement  of  vanes 
consisted  of  blades  made  from  one  quarter  cross  sections  of 
1.625  inch  O.D.  copper  tubing,  with  a  wall  thickness  of  0.04 
inches .  Five  blades  were  used  and  they  were  spaced  so  that 
the  distance  between  them  increased  in  arithmetic  progression, 
starting  at  the  upper  corner  as  shown  in  Figure  3. 

To  solve  the  problem  of  obtaining  symmetrical  flow 
at  the  tube  inlet,  the  one-inch  diameter  pipe  between  the  tee 
and  the  entrance  cone  was  increased  to  40  pipe  diameters 
length.  At  this  length,  within  the  one-inch  tube  the  flow  of 
its  own  accord  could  be  expected  to  become  fully  developed, 
and  thus  to  produce  a  symmetrical  velocity  profile. 

At  the  same  time,  it  was  found  desirable  to  perfect 
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the  connection  from  the  entrance  cone  to  the  six- inch  tube. 

This  new  connection  is  shown  in  Figure  4. 

The  final  form  of  the  entrance  section,  shown  in 
Figure  5,  gave  symmetrical  velocity  profiles  at  all  points 
along  the  axis  of  the  aluminum  tube  as  verified  by  the 
experimental  data, 
b .  The  Air  Supply 

Air  was  supplied  by  a  Gardner-Denver  Model  3CDL-5 
Cyclo  blower,  a  rotary  lobe  type  of  axial  flow  blower  pro¬ 
ducing  a  continuous  pulsat ion-f ree  discharge.  The  blower 
was  direct  driven  by  a  20  hp,  3600  rpm  induction  motor,  which 
had  a  maximum  capacity  of  240  SCFM  at  12  psig.  The  compression 
chamber  was  so  designed  that  lubrication  was  unnecessary,  thus 
guaranteeing  an  oil-free  air  supply.  The  operation  and 
maintenance  of  this  unit  is  described  in  the  Appendix  section. 

The  blower  intake  was  protected  by  an  air  filter. 

A  relief  valve  was  installed  in  the  discharge  line  to  pre¬ 
vent  possible  over  pressure.  One  check  valve  was  installed 
to  protect  the  filter  and  another  to  prevent  back-flow  of 
bulk  material  into  the  blower  at  shutdown. 

Because  of  the  blower  air  discharge  temperatures, 
between  200  and  300°F,  an  aftercooler  was  required.  The 
final  air  temperatures  were  measured  with  an  iron  constantan 
thermocouple  connected  to  a  Brown  single-point  temperature 
recorder  .  By  suitable  manual  control  of  the  cooling  water  flow 
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rate,  the  aftercooler  not  only  cooled  the  air  but  also  main¬ 
tained  it  at  a  constant  reference  temperature  for  all  experi¬ 
mental  runs. 

The  temperature  recorder  was  checked  with  a  calibrated 
mercury  thermometer.  At  operating  temperatures,  the  deviation 
was  less  than  0.3°F,  and  at  steam  and  freezing  point  slightly 
less  than  1°F.  Due  to  the  fact  that  the  cooling  water  had  to 

be  controlled  manually,  the  air  temperature  could  only  be 

.  .  + 

guaranteed  within  -1°F.  A  startup  period  of  about  one  hour 
was  required  to  obtain  a  constant  air  temperature. 

A  moisture  separator  was  installed  downstream  from 
the  aftercooler  to  eliminate  entrained  water. 

The  air  flow-rate  was  controlled  by  an  orifice  flow¬ 
meter  of  diameter  of  1.544  inches,  which  gave  sufficient  control 
over  the  entire  flow  range.  Two  u— tube  manometers  were  con¬ 
nected  across  the  orifice,  one  with  water  as  the  manometer 
fluid  for  low  pressure  differentials  and  the  other  with  mer¬ 
cury  for  higher  differentials.  Both  manometers  could  be  used 
for  measurement  of  pressure  differentials  across  the  orifice 
or  for  gage  pressure  on  the  upstream  side  of  the  orifice. 

The  flow  rate  was  set  manually  by  means  of  a 
pneumatically  operated  1-1/4  inch  parabolic  needle  valve, 

Foxboro  type  F-8 ,  located  in  the  2-inch  recycle  line.  To  ob¬ 
tain  very  low  flow  rates,  it  was  necessary  to  open  the  globe 
valve  in  the  bleed-off  line. 
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In  measuring  air  flow  with  the  orifice  meter,  the 
geometry  of  the  2-1/2  inch  air  line  from  the  blower  to  the 
entrance  cone  precluded  having  an  adequate  straight  run  up¬ 
stream  from  the  orifice.  To  eliminate  the  problem  of  a 
fluctuating  orifice  differential  pressure  which  resulted, 
straightening  vanes  were  designed  and  installed  48  inches  up¬ 
stream  according  to  American  Gas  Association  specif ications ( 2 3) . 
The  fluctuations  were  thus  reduced  to  negligible  values  for 
all  flow  rates  encountered, 
c .  The  Velocity  Profile  Probe  Mount 

A  probe  mount  arrangement  was  available! 1)  to  permit 
the  positioning  of  probes  along  the  tube  diameter  at  right 
angles  to  the  axis  of  the  tube.  The  position  of  the  probe 
could  be  read  from  a  vernier  scale  with  0.01  inch  subdivisions. 
At  some  of  the  probe  locations,  the  orientation  of  the  probe 
stem  travel  could  be  varied  by  90  degrees.  This  allowed  the 
measurement  of  two  velocity  profiles  at  right  angles  to  each 
other . 

Because  a  total-head  pitot  tube  was  employed  to 
measure  mean  point  velocity,  the  absolute  pressure  at  the 
probe  location  and  the  difference  between  total  pressure  and 
static  pressure  at  the  same  location  were  required.  The  pitot 
tube  consisted  of  a  23-gauge  hypodermic  needle  as  shown  in 
Figure  6.  The  static  pressure  was  measured  at  the  wall  of 
the  tube.  A  very  fine  hole,  0.031  inch  diameter,  was  drilled 
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in  a  brass  plug  that  had  been  machined  to  fit  with  the  inner 

surface  of  the  tube.  This  tap  arrangement  is  shown  in  Figure 

«  * 

7 .  The  static  gauge  pressure  was  read  on  a  u-tube  manometer 
using  water  as  the  manometer  fluid.  Pressure  differentials 
were  measured  with  a  Flow  Corporation  Model  MM3  Micromanometer. 
The  micromanometer  consisted  of  a  u-tube  type  manometer  operat¬ 
ing  on  a  nulling  principle.  It  measured  a  maximum  pressure 
difference  of  two  inches  of  manometer  fluid  with  a  resolution 
of  io.OOOl  inch  and  an  accuracy  of  -0.0002  inch  of  manometer 
fluid.  In  the  present  work,  the  manometer  fluid  employed  was 
n-Butyl  alcohol  with  a  specific  gravity  of  0.8166  at  75°F. 

The  zero  setting  changed  somewhat  with  time,  and  it 
was  found  necessary  to  measure  and  record  the  zero  setting 
between  each  pressure  differential  measured  and  use  the  average 
value . 

The  barometric  pressure  was  obtained  with  a  standard 
barometer  and  was  corrected  for  variations  in  the  ambient 
temperature,  an  estimated  accuracy  of  io.2  mm  Hg  was  obtained, 
d .  The  Center-Line  Velocity  Travelling  Carriage 

A  carriage  arrangement  with  pitot  tube  and  static 
pressure  probe  was  designed  to  obtain  measurements  of  the 
center-line  velocity  at  any  axial  point  along  the  entrance 
section  of  the  aluminum  tube.  The  carriage,  as  shown  in 
Figure  9,  consisted  of  a  short  piece  of  7/8  inch  pipe  with 
six  arms,  mounted  three  on  each  end,  to  secure  the  central 
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positioning  of  the  carriage.  The  tip  of  each  arm  was  fatted 
with  a  small  teflon  bearing  surface  to  protect  the  surface 
of  the  test  pipe  and  to  make  the  carriage  easier  to  move.  To 
obtain  the  total  head  at  the  center-line  of  the  tube,  a  hypo¬ 
dermic  needle  was  mounted  on  the  front  end  of  the  7/8  inch  pipe 
at  the  center  of  the  carriage.  A  static  pressure  needle  was 
mounted  on  one  of  the  arms,  one-third  tube  radius  distant  from 
the  tube  wall. 

The  carriage  could  be  positioned  along  the  tube  axis 
by  a  steel  cable  connected  to  the  carriage  and  which  extended 
up  and  out  of  the  tube  through  the  upper  cone.  The  leads  from 
the  total  head  and  the  static  pressure  probes  consisted  of  0.1 
inch  O.D,  polyethylene  tubing  taped  to  the  steel  cable  to 
ensure  rigidity.  The  arrangement  is  shown  in  Figure  10. 
e •  The  Measurement  of  Pressure  Drop  Along  the  Tube 

Pressure  drops  along  the  tube  were  measured  with  two 
wall  static  pressure  taps  of  the  type  shown  in  Figure  7 . 
Positioned  a  known  distance  apart,  the  taps  provided  a  means 
for  obtaining  pressure  drop  data  between  any  two  of  the  fixed 
probe  mount  locations. 
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2 .  Experimental  Procedure s 
a •  Calibration  of  the  Orifice 

The  orifice  was  calibrated  by  integrating  velocity 
profiles  obtained  with  the  pitot  tube  as  the  standard.  The 
mean  point  velocities  across  the  tube  were  measured  and  then 
integrated  to  give  an  average  or  bulk  velocity, 

1  r 

uav  “  — ~2  /  u  ^irrdr  (51) 

TTR  J 

The  flow  rate  in  standard  cubic  feet  per  minute,  Qq,  was  ob¬ 
tained  from  the  following  equation, 

T  p  1 

Qo  =  uav  F  —  — ‘  —  <52) 

T  PQ  60 

When  the  temperature  of  the  air  supply  is  held  constant  at 
531. 6°R,  equation  (52)  may  be  simplified  to 


Qq  =  0.0167  Puav  (53) 

The  flow  rate  can  also  be  expressed  in  terms  of  the  pressure 
drop  across  the  orifice 

A  I  2g  AP 

Q  =  CFYM  — —  (54) 


Upon  substitution  of  the  cross-sectional  area  of  the  tube  and 
converting  the  units  to  standard  cubic  feet  of  air  per  minute, 
one  obtains 


AP  P 


(  55) 


31,73  CY 


T 
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An  average  value  of  the  orifice  coefficient  was  obtained  by 
plotting  Qq  from  equation  (53)  versus  Qq/C  from  equation  (55) 
and  measuring  the  slope  of  the  resulting  straight  line.  Except 
for  small  deviations,  the  calibration  points  fell  on  a  straight 
line  passing  through  the  origin.  The  flow  coefficient,  obtained 
by  the  method  of  least  squares,  was  equal  to  0.64  with  a  maxi¬ 
mum  deviation  of  3%.  The  correlation  is  shown  in  Figure  11, 
and  experimental  data  are  given  in  Table  3.  The  expansion 
factors,  Y,  which  were  applied  ranged  from  0.97  to  1.0. 
b .  Measurement  of  Velocity  Profiles 

Complete  velocity  profiles  were  measured  at  a  dis¬ 
tance  of  41  pipe  diameters  from  the  entrance,  over  a  Reynolds 
number  range  of  9,000  to  56,000,  using  the  pitot  tube  together 
with  the  wall  static  pressure  tap.  Values  of  the  mean  point 
velocities  were  obtained  at  positions  as  close  as  0.0775  inch 
to  the  tube  wall. 

At  Reynolds  numbers  less  than  9,000,  the  impact 
pressure  was  too  small  to  be  measured  accurately  by  the  micro¬ 
manometer,  The  upper  limit  of  the  Reynolds  number,  56,000, 
corresponded  to  the  capacity  output  of  the  blower. 

Two  traverses  were  taken  at  right  angles  to  each 
other  at  distances  of  41  and  also  4,2  diameters,  downstream 
from  the  entrance,  to  ensure  that  the  velocity  profiles  were 
symmetrical  about  the  tube  axis.  These  profiles  are  shown  for 
Re  =  18,000  at  41  diameters,  and  for  Re  =  24,000  at  4.2  diameters, 
downstream  from  the  entrance  in  Figure  12.  The  corresponding 
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experimental  data  are  given  in  Tables  5  and  6. 

The  velocity  profiles  corresponding  to  Re  =  13,000 
measured  at  35.2  and  repeated  at  41  diameters  were  identical 
within  the  limitations  of  the  experimental  measurements,  hence, 
it  was  concluded  that  the  flow  was  fully  developed  at  the  41 
diameter  position. 

c .  Measurement  of  Center-Line  Velocity 

To  obtain  a  profile  of  axial  velocity  development 
using  the  travelling  carriage  already  described,  certain  dif¬ 
ficulties  may  be  anticipated.  The  presence  of  the  carriage 
could  be  expected  to  influence  the  static  pressure  in  the 
tube.  The  static  pressure  was  measured  about  four  inches 
ahead  of  the  main  structure  of  the  carriage.  Initial  tests 
comparing  the  static  pressure  obtained  with  the  carriage  static 
pressure  probe  and  that  of  a  wall  tap,  having  the  carriage 
a  few  inches  further  downstream  from  the  location  of  the  tap, 
showed  no  influence  of  the  carriage.  However,  the  velocity 
obtained  with  the  carriage  was  consistently  higher  than  that 
obtained  with  the  traversing  pitot  tube  when  both  were  measured 
at  the  same  location.  Assuming  that  differences  between  the 
total  pressures  are  minor  conpared  to  those  between  the  static 
pressures  and  that  the  carriage  values  contain  the  error,  it 
could  be  deduced  that  the  static  pressure  obtained  with  the 
carriage  read  too  low. 

Experiments  were  then  undertaken  using  the  static 
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pressure  taps  along  the  tube  wall  together  with  the  carriage 
pitot  tube.  The  carriage  was  so  positioned  that  the  static 
pressure  hole  at  the  wall  was  two  inches  upstream  from  the 
center-line  pitot  tube.  Center-line  velocities  measured  in 
this  way  were  noticeably  lower,  but  these  values  were  still 
higher  than  those  obtained  in  the  corresponding  velocity  pro¬ 
files.  It  was  thus  concluded  that  the  static  pressure  was 
affected  by  the  presence  of  the  carriage. 

The  results  of  carriage-measured  center-line 
velocities  obtained  by  using  both  the  static  pressure  probe  at 
the  carriage  and  the  static  pressure  taps  at  the  tube  wall 
are  compared  in  Figure  13.  The  data  are  given  in  Tables  38  to 
43. 

d ,  Measurement  of  Friction  _Fa ct or s 

The  static  pressure  probes  shown  in  Figure  8,  spaced 
five  feet  apart  at  44.8  and  35.2  diameters  downstream  from  the 
entrance,  were  used  to  measure  the  pressure  drop  along  the 
tube.  Blasius  friction  factor  were  calculated  assuming  in¬ 
compressible  flow, 

2g  D  AP 

f  =  -----  (-  — )  (56) 

2  Ax 

P  uav 

The  friction  factors  thus  obtained  disagreed  greatly  with  the 
expected  values (6) ,  Furthermore,  by  interchanging  the  two 
probes,  completely  different  values  were  obtained.  It  was 
then  decided  to  use  wall  taps  to  measure  pressure  drop  for 


, 
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the  calculation  of  friction  factor.  Although  the  wall  taps 
were  no  more  mechanically  perfect,  this  arrangement  was  believed 
to  be  less  sensitive  to  mechanical  imperfection  since  the  air 
velocity  close  to  the  wall  was  low. 

The  friction  factors  so  obtained,  given  in  Table  4, 
are  plotted  and  compared  to  the  Blasius  friction  factor  cor¬ 
relation  in  Figure  14.  The  Blasius  correlation,  may  be  considered 
for  all  practical  purposes  to  be  identical  to  that  of  Moody(6), 
in  the  Reynolds  number  range  investigated.  The  Blasius  friction 
factor  correlation  is  represented  by 


f  = 


0.316  Re 


-1/4 


(57) 


The  data  so  obtained  were  slightly  higher  than  those  given  by 
equation  (57).  The  type  of  deviation  encountered  did  not 
seem  to  be  caused  by  ordinary  roughness  but  was  more  of  the 
type  caused  by  wall  waviness.  Prandtl  and  Tietjens(4)  dis¬ 
cuss  the  two  types  of  roughness  and  also  mention  that  a  com¬ 
bination  of  roughness  and  waviness  has  been  encountered  with 
drawn  metal  tubes  such  as  the  aluminum  tube  involved  in  these 


tests . 

The  outside  diameter  of  the  tube  was  measured  and 
found  to  vary  slightly.  The  maximum  variation,  3/32  inch,  was 
found  in  the  upper  region  of  the  test  tube  at  approximately 
45  diameter  from  the  entrance  while  along  the  remainder  of  the 
tube,  variations  of  less  than  1/32  inch  were  encountered. 
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The  data  so  obtained  were  fitted  by  the  method  of 
least  squares  and,  assuming  waviness  to  be  present,  are  re¬ 
presented  by 

f  =  0.377  Re“1//4  (58) 

However,  since  the  maximum  variation  in  diameter  occurred 
downstream  from  the  location  where  the  fully  developed  velocity 
profiles  were  obtained,  the  Blasius  friction  factor  correlation 
was  used. 
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IV.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


1.  Fully  Developed  Flow 


Velocity  profiles  were  obtained  at  41  diameters 


from  the  entrance  (measured  from  the  junction  of  the  inlet  cone 
and  the  6.255  inches  tube)  for  flow  rates  corresponding  to  8 
different  Reynolds  numbers.  Values  of  the  mean  point  velocity 
were  calculated  from  the  following  equation: 


(59) 


u 


No  pitot  tube  calibration  constant  was  applied  in  equation  (59) 
since  the  static  pressure  was  measured  at  the  tube  wall . 
Measured  and  calculated  point  velocity  data  are  given  in  Tables 
7  to  14.  For  all  profiles  measured,  the  point  velocity  was 
obtained  on  both  sides  of  the  center-line,  the  two  values 
being  distinguished  in  the  Tables  by  "in"  and  "out".  The 
value  of  the  velocity  used  in  the  correlations  was  the  average 
of  "in"  and  "out".  The  velocity,  u,  is  plotted  versus  radial 
position,  r,  in  Figure  15.  To  obtain  an  average  velocity, 
equation  (51)  was  integrated  by  means  of  a  planimeter .  The 
ratios  of  average  to  maximum  velocity  was  computed  and  plotted 
versus  Reynolds  number  on  Figure  16. 
a.  The  Power  Law 


The  point  velocities  were  correlated  in  terms  of 


the  power  law  given  by  equation  (1)  and  were  found  to  correlate 
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well  if  the  points  measured  close  to  the  wall,  at  y  =  0.0775, 
were  omitted.  The  next  value  was  obtained  at  y  =  0.1275,  and 
was  found  to  be  more  reliable.  The  data  given  in  Tables  15 
to  18  correlated  with  an  average  deviation  of  1.8%  and  of  3.4% 
when  the  point  close  to  the  wall  was  included.  The  power  law 
exponents  so  obtained  are  shown  in  Table  19,  and  are  plotted 
versus  Reynolds  number  in  Figure  17.  Nunner(19)  reported  that 
the  power  law  exponent  could  be  expressed  in  terms  of  the 
friction  factor 


n 


(60) 


Upon  combining  this  with  the  Blasius  friction  factor  of 
equation  (57),  the  exponent  can  be  related  to  the  Reynolds 
number,  giving 


n 


(61) 


Equation  (61)  correlates  the  experimental  data,  as  shown  in 
Figure  17,  reasonably  well. 

4-  *l-i 

The  dimensionless  velocity  raised  to  the  n  power, 


u  n  y 

( — ) ,  was  plotted  versus  —  ,  using  the  value  of  n  from  equation 


u 


m 


R 


(61) .  It  was  found  that  the  data  fell  on  a  straight  line 


(~“)n 
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R 


(62) 
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except  for  the  center-line  velocity.  The  average  value  of  C 
was  found  by  least  squares  to  equal  1.145  and  it  varied  by 
less  than  2.2%.  The  results  are  shown  in  Figure  18.  To  avoid 
crowding  of  the  data  on  the  graph,  the  ordinates  have  been 
shifted  for  each  run. 

By  integrating  the  Power  Law,  Schlichting(6)  obtained 

an  expression  for  the  ratio  of  average  to  maximum  velocity, 

.2 

- -  (63) 


u 


av 


2n‘ 


u 


m 


(n+1) (2n+l) 


Both  equation  (63)  combined  with  equation  (61),  and  equation 
(43)  gave  a  velocity  ratio  slightly  lower  than  the  experimental 
values.  The  experimental  results  of  Rothfus  and  Monrad(8), 
however,  agree  closely  with  those  obtained  in  this  work.  The 
results  given  in  Table  19  are  shown  and  compared  in  Figure  16. 
b .  The  Logarithmic  Velocity  Distribution 

In  the  preceding  theory,  equations  (8)  and  (12) 
were  left  with  one  constant,  K,  to  be  determined  experimentally. 
It  was  also  assumed  to  be  a  universal  constant,  that  is  to  say, 
not  dependent  on  Reynolds  number,  A  value  of  I<  =  0.40  has 
been  experimentally  obtained  by  Nikuradse( 6) ,  but  with  the 
present  data,  considerably  lower  constants  were  obtained.  By 
the  method  of  least  squares  applied  to  all  points  measured, 
the  following  values  of  K  were  obtained: 

K  =  0.348  ,  for  equation  (8) 

K  =  0.311  ,  for  equation  (12) 


0-t^ 

' 


For  the  individual  runs,  K  varied  between  0.28  and  0.38  with 
no  consistent  decrease  or  increase  with  increasing  Reynolds 
number  and  with  average  deviations  of  more  than  10%.  The 
experimental  data  and  calculated  values  are  given  in  Tables 
20  to  23.  To  avoid  crowding  the  experimental  results  for 
Reynolds  numbers  of  56,045  and  18,176  only,  together  with 
equations  (8)  and  (12),  using  the  two  experimental  values  of 
K  and  the  value  of  0.4  are  shown  in  Figure  19.  On  the  basis 
of  the  data  obtained  in  the  present  investigation,  the  cor¬ 
relation  of  Prandtl,  equation  (8),  was  superior  to  that  of 
von  Karman,  equation  (12),  but  both  were  inferior  to  the  simple 
Power  Law  distribution. 

An  equation  similar  to  the  von  Karman  equation  can 
be  obtained  if  we  neglect  to  apply  the  second  boundary  con¬ 
dition,  u  =  um  for  y  =  R.  Using  this  equation,  in  which  there 
are  two  constants  to  be  evaluated,  a  considerably  better  fit 
was  achieved.  However,  one  set  of  constants  was  needed  for 
each  value  of  Reynolds  number . 

The  correlation  proposed  by  Wang,  equation  (13), 
was  not  evaluated  since  it  was  too  inconvenient  for  practical 
use . 

c.  The  Universal  Velocity  Distribution 

Equation  (20)  provides  the  most  widely  used  basis 
for  correlating  point  velocities  in  a  duct,  at  a  distance  from 
the  wall.  The  equation  contains  two  constants  which  were 
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determined  experimentally  and  found  to  be  K  =  0.361  and  C'  =  4.00 
by  the  method  of  least  squares.  Equation  (20),  as  obtained  by 


Nikuradse,  Deissler 

and 

in  the  present  work  are, 

respectively, 

u+ 

=  2 

.5  In  y+  +  5.5 

(23) 

u+ 

=  2 

.78  In  y+  +  3.8 

(24) 

U+ 

=  2 

.77  In  y+  +  4.0 

(64) 

Based  on  Nikuradse’ s  data,  equation  (23)  was  found  to  be  valid 

JL 

for  y'  >30,  while  Deissler  reported  equation  (24)  to  be  valid 
for  y  >26.  The  present  work  seems  to  agree  with  Deissler. 

The  great  difference  in  constants  given  that  of  equation  (23) 
and  those  of  (24)  and  (64)  is  believed  to  have  been  caused  by 
the  different  range  of  Reynolds  numbers  over  which  they  were 
obtained.  Nikuradse  worked  with  Reynolds  numbers  ranging  from 
4,000  to  3,240,000  while  Deissler,  and  to  an  even  greater  ex¬ 
tent  the  present  investigator,  worked  over  a  considerably 
smaller  range  of  Reynolds  numbers.  It  is  known  that  all  of 
the  points  of  a  complete  velocity  profile  do  not  fall  on  a 
straight  line  as  given  by  equation  (20),  a  true  curve  giving 
slightly  higher  velocities  near  the  center  line.  This  effect 
is  confirmed  on  Figure  20.  The  present  work  and  Deissler' s 
did  not  include  velocity  profiles  in  the  range  of  Reynolds 
number  from  4,000  to  9,000,  and  thus  the  effect  of  higher 
point  velocities  than  the  equation  (20)  indicates  was  not 
present  in  the  lower  region,  y+<200.  This  effect  could  account 
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for  a  steeper  slope  and  a  shorter  intersection  obtained  in 
equations  (23)  and  (64).  Equations  (23)  and  (64)  together 
with  the  experimental  results  are  shown  in  Figure  20  and  the 
corresponding  data  are  given  in  Tables  24  to  27.  Deissler's 
correlation  was  omitted  from  Figure  20  because  it  was  prac¬ 
tically  identical  to  that  of  equation  (24).  For  y+  >  26,  an 
average  percent  deviation  of  1.28%  and  a  maximum  deviation  of 
of  3.96%  were  encountered  for  the  62  points  which  were  measured. 
Equation  (64)  with  one  set  of  constants  gave  a  better  fit  over 
the  entire  range  of  Reynolds  numbers  than  the  von  Karman 
equation  using  a  different  set  of  constants  for  each  in¬ 
dividual  Reynolds  number . 

By  examining  the  constants  obtained  for  equation  (20) 
in  the  individual  runs,  it  was  found  that  they  could  vary  greatly 
without  causing  much  change  in  the  computed  velocity.  Further¬ 
more,  it  was  noticed  that  for  a  decrease  in  slope  the  constant 
representing  the  intersection  always  increased  and  vice-versa. 

The  maximum  variation  of  K  for  any  individual  run  was  10% 
while  C*  had  a  maximum  deviation  of  more  than  30%.  This  sur¬ 
prisingly  high  variation  in  the  value  of  the  constants  only 
corresponded  to  a  maximum  variation  of  4%  in  computed  velocity. 
Errors  in  the  constants  attributed  to  possible  errors  in 
the  experimental  measurements  giving  the  velocity  profile  will 
be  analyzed  in  the  appendix  section.  The  results  of  the 
analysis  clearly  show  that  the  value  of  the  constant  C'  was 
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far  more  sensitive  than  the  value  of  K  to  errors  in  the 
variables  u'  and  y+ . 

To  account  for  the  observed  deviations  from  equation 
(64),  an  equation  of  the  form  suggested  by  Bogue,  equation  (35) 
was  assumed 

+  Y  + 

u  -  C  (-- , f)  =  A  In  y  +  B  (65) 

R 


where 


C  (— ,  f) 
R 


exp  ( 


(66) 


The  constants  A  and  B  were  obtained  by  the  method  of  least 
squares  while  the  constants  C,  D  and  E  were  determined  by  a 
trial-and-error  procedure  as  shown  in  the  Appendix.  The 
criterion  used  for  a  best  fit  was  the  minimum  of  the  sum  of 
the  absolute  values  of  percent  error .  The  constants  which 
were  obtained  were  slightly  different  from  those  of  Bogue  who 
selected  the  constants  A  and  B  in  equation  (65)  so  that  this 
equation  would,  upon  integration,  yield  the  Moody  friction 
factor  correlation  for  smooth  pipes  and  then  curve  fitted  the 
remaining  constants.  The  present  work  gave 

u+  -  C  (-,  f)  =  2.43  In  y+  +  5.24  (67) 
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On  the  basis  of  previous  findings,  equation  (67)  was  assumed 
to  be  valid  for  y'  26.  An  average  percent  deviation  of 
0.81%  and  a  maximum  deviation  of  4.5%  resulted. 

However,  by  examining  the  data  obtained  in  individual 
runs,  it  was  found  that  the  effect  of  Reynolds  number  on  the 
velocity  profiles  was  nearly  insignificant,  if  at  all  observ¬ 
able.  Thus,  the  empirical  correction  function,  equation  (68), 
was  simplified  as  follows: 

_(2  _  D')2 

C  (-)  =  C'  exp  ( - 2 - )  (69) 

R 

E  ' 

Using  the  same  procedure  as  previously,  the  following 
correlation  was  obtained: 
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(71) 


The  correction  function  given  by  equation  (71)  is  of  simpler 
form  than  equation  (68)  because  it  can  be  represented  in  the 
form  of  a  single  curve,  a  function  of  radial  position  y/R . 
Over  the  limited  range  investigated,  Reynolds  number  from 
9000  to  56,000,  equation  (70)  v/as  found  to  be  comparable  to 
equation  (67) .  An  average  percent  deviation  of  0.82%  and  a 
maximum  deviation  of  4.26%  were  encountered.  The  correction 
function,  equation  (71),  is  shown  in  Figure  21.  The  experi- 
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mental  data  are  given  in  Tables  24  to  27  and  equation  (70)  is 
plotted  in  Figure  22 . 

d .  The  Velocity  Distribution  in  the  Vicinity  of  the  Wall 

The  limited  experimental  data  which  were  obtained  in 
the  region,  y 1  <  26,  showed  considerable  scatter  as  may  be 
seen  in  Figure  20.  With  only  9  points,  and  a  maximum  of  2 
points  obtained  in  a  given  run,  no  conclusion  may  be  drawn  as 
to  whether  one  correlation  fits  better  than  another.  However, 
the  points  indicate  that  there  is  a  changed  physical  situation 
for  y+  <  25.  This  change  is  caused  by  the  wall  which  super¬ 
imposes  a  viscous  shear  stress  on  the  fluid  near  the  wall. 

The  experimental  data  reported  here  are  insufficient  to  deter¬ 
mine  the  exact  limit  where  the  viscous  shear  stress  effect  on 
the  flow  diminishes.  To  describe  this  situation  more  vividly, 
the  empirical  correlation,  equation  (22),  together  with  that 
of  Deissler,  equation  (31),  and  that  of  Wasan  et  al,  equation 
(33),  are  shown  for  comparison  in  Figure  20. 

Velocities  computed  from  equation  (22)  together 
with  (23),  and  from  (31)  together  with  (64),  will  not  give 
continuous  velocity  profiles.  Equation  (33)  was  so  derived 
that  the  resulting  velocity  profile  would  have  a  continuous 
transition  to  equation  (23),  however,  it  is  unlikely  that 
equation  (23)  is  valid  down  to  y+  =20.  By  using  the  approach 
of  Wasan  et  al(12),  matching  equation  ( 32)  with  equation  (20), 
but  using  the  constants  obtained  in  the  present  work,  the 
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following  relation  was  found: 

u+  =  y+  -  1.403  (10-4)(y+)4  +  4.55  (10-°)(yVJ  (65) 

The  transition  was  found  to  occur  at  y+  =  18.  As  can  be  seen 
in  Figure  20,  the  above  equation  agrees  more  closely  with  the 
equation  (22),  that  has  been  experimentally  verified  in  the 
region.  The  derivation  of  equation  (65)  is  given  in  the 
Appendix  section.  The  computed  values  are  given  in  Table  28. 
Equation  (65)  possesses  the  advantages  of  providing  a  continuous 
transition  to  the  equation  valid  in  the  turbulent  core,  of 
approaching  laminar  flow  for  y+  =  5,  and  it  is  more  convenient 
to  use  than  the  Deissler  correlation,  equation  (31).  However, 
experimental  confirmation  is  lacking. 

Gill  and  Scher  '  s  velocity  distribution  equation  (39), 
based  on  a  modified  Prandtl  mixing  length  theory  was  tested 
with  the  present  data.  The  Legendre-Gauss  quadrature  formula (20) 
was  used  for  the  numerical  integration.  Details  of  the  com¬ 
putational  procedure  are  given  in  the  Appendix  section.  The 
velocity  computed  from  equation  (39)  resulted  in  a  curved 
velocity  profile  with  semilogarithmic  coordinates  as  can  be 
seen  in  Figure  23.  The  distribution  which  was  obtained  was 
consistently  too  low  in  the  region  close  to  the  center  line. 

The  two  constants  which  were  employed,  determined  by  Gill  and 
Scher,  were  varied  but  no  significant  improvements  were  achieved. 
For  computations,  y  was  obtained  from  the  following  relation: 
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The  computed  values  are  given  in  Tables  29  to  32 .  The  cor¬ 
relation  follows  equation  (21)  and  (22)  closely  in  the  region 
near  the  wall .  In  the  region  far  away  from  the  center  line  it 
is  in  fair  agreement  with  experimental  results  and  equation  (64), 
however,  approaching  the  center  line  the  correlation  is  not  very 
satisfactory.  Computed  values  are  given  in  the  Appendix  section 
following  the  computational  procedure.  Experimental  data  and 
computed  values  are  plotted  in  Figure  23  for  Reynolds  numbers 
of  13,120  and  56,050. 


The  Pai  power  series  correlation,  given  by  equation 


(42),  was  found  to  be  superior  to  the  other  correlations  valid 
from  the  v/all  to  the  center  line.  By  means  of  the  Brodkey 
correlation ( 15)  for  the  variables  used  in  the  Pai  type  equation, 
S  and  M,  the  velocity  profiles  shown  in  Figure  24  were  obtained. 
The  large  deviations  for  points  measured  close  to  the  wall 
were  random  and  thus  believed  to  have  been  caused  primarily  by 
inaccuracies  in  measuring  the  velocity  and  not  by  an  incon¬ 
sistency  in  the  correlation.  Points  which  were  obtained  in 
the  turbulent  core,  y+  >  26,  showed  an  average  deviation  of 
only  0.9%  and  a  maximum  deviation  of  4.5%.  This  was  consider¬ 
ably  better  than  what  was  obtained  with  equation  (64) ,  and 
comparable  to  the  results  obtained  with  the  Bogue  universal 
velocity  distribution,  equation  (35).  Furthermore,  for  the 
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Pai  power  series  correlation  no  constants  were  curve-fitted 
to  the  present  experimental  data.  Brodkey(15)  reported  certain 
discrepancies  at  high  Reynolds  number.  In  the  range  tested, 
9,000  <  Re  <  56,000,  no  such  trend  was  observed.  For  decreas¬ 
ing  Reynolds  numbers,  both  the  experimental  data  and  equation 
(42)  converge  to  the  laminar  velocity  distribution,  as  can  be 
seen  in  Figure  24. 

The  von  Karman  correlation  was  re-arranged  in  terms 
of  u/um  versus  y/R  and  plotted  in  Figure  24.  Dividing  equation 
(12)  by  um  gave 


where 


and 


1 

1  +  — 


K 


u* 


In  (1 


um  (n+1)  (2n+l) 
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Equations  (68)  and  (63)  were  substituted  into  equation  (67) 
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The  power  lav/  exponent,  n,  is  related  to  the  Reynolds  number 
via  equation  (61),  the  Blasius  friction  factor  via  equation  (57) 
and  finally  the  constant  K  was  determined  for  each  individual 
run  by  the  method  of  least  squares.  The  correlation  represented 
the  velocity  profiles  very  well,  and  provides  an  interesting 
comparison  to  the  Pai  correlation.  At  low  Reynolds  number, 
the  curvatures  seemed  to  approach  each  other.  At  higher  Rey¬ 
nolds  number  they  differ  considerably  for  y/R  between  0.2  and 
0.5  and  the  experimental  data  tend  to  fall  between  the  tv/o 
curves.  Good  agreement  between  experimental  results  and 
equation  (69)  could  only  be  obtained  by  using  a  different  I< 
for  each  individual  run.  The  K-value  varied  by  more  than  10% 
with  varying  Reynolds  number,  but  no  consistency  could  be 
observed . 

This  type  of  a  correlation,  equation  (69),  could 
only  be  useful  if  K  proved  to  be  a  universal  constant,  indepen¬ 
dent  of  Reynolds  number,  and  the  value  of  K  could  be  deter¬ 
mined  accurately.  In  the  present  work  the  average  value  of 
K  was  0.31.  Data  calculated  for  the  Pai  power  series  and  for 
equation  (69)  are  given  in  the  Appendix  section  in  Tables 
3  3  to  36 . 


2  . _ The  Entrance  Section 

Velocity  profiles  were  obtained  at  4.2  diameters 
downstream  from  the  entrance  for  flow  rates  corresponding  to 
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those  measured  at  41  diameters.  The  entrance  velocity  profiles 
are  shown  in  Figure  25,  and  the  corresponding  values  are  given 
in  Table  37 .  The  dimensionless  velocity,  u/um,  was  plotted 
versus  y/R  in  Figure  26  for  Reynolds  number  18,000,  29,000  and 
38,000  and  both  4.2  and  41  diameter  locations.  Upon  examining 
the  two  sets  of  profiles  shown  in  Figure  26,  it  can  be  seen 
that  the  velocity  profiles  obtained  near  the  entrance  are 
fuller,  i.e.  flatter  near  the  center  line,  than  those  obtained 
at  41  and  also  that  both  sets  of  velocity  profiles  become 
fuller  with  increasing  Reynolds  number.  The  velocity  profiles 
in  Figure  26  represent  the  limiting  conditions  for  the  test 
unit  and  the  evaluation  of  the  change  from  the  one  to  the 
other  was  one  of  the  aims  in  this  work. 

To  provide  a  comparison  of  the  velocity  distributions, 
the  entrance  profiles  were  tested  with  the  universal  velocity 
distribution  correlation,  equation  (20) .  The  constants,  K  and 
C ' ,  were  obtained  by  the  method  of  least  squares  giving  the 
equation 

u+  =  1.68  In  y+  +  9.77  (70) 

This  correlation  together  with  the  experimental  points  are 
shown  in  Figure  27 .  All  of  the  points  measured  were  included 
and  except  for  the  data  obtained  at  the  lowest  flow  rate, 

Re  =  8,900,  the  data  correlated  reasonably  well  with  an 
average  deviation  of  1.39%  and  a  maximum  deviation  of  4.44%. 
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The  question  arises  whether  a  set  of  profiles  at  any 
fixed  location  with  varying  Reynolds  number,  could  be  plotted 
on  a  single  curve.  Since  the  universal  velocity  distribution 
makes  no  distinction  between  different  fully  developed  profiles 
in  turbulent  flow,  this  approach  was  attempted,  however,  the 
question  could  not  be  solved  with  the  limited  data  available. 

Equation  (70)  intersects  the  correlation  for  fully 

developed  flow  at  y1"  =  2  00.  The  slope  of  the  equation,  1.68, 

indicates  much  flatter  profiles  than  for  fully  developed  flow. 

For  a  pipe  with  an  entrance  giving  constant  velocity  across 

the  diameter  (plug  flow) ,  a  velocity  profile  would  yield  a 

horizontal  line  in  terms  of  the  universal  velocity  distribution 

correlation.  If  the  velocity  profiles  at  various  locations 

between  the  entrance  and  the  location  where  the  flow  becomes 

fully  developed  would  correlate  in  terms  of  straight  lines 

on  a  u*  -  y+  plot,  their  slopes  should  vary  between  2.77  and 

zero.  In  this  way,  the  universal  velocity  distribution  could 

be  used  for  correlating  data  in  the  entrance  section  of  a 

pipe.  However,  it  is  apparent  that  a  correction  factor  of 

the  type  suggested  by  Rothfus  et  al(8,9)  would  be  necessary 

in  a  "universal"  correlation  because  the  intersection  for 

flat  profiles  would  increase  with  increasing  Reynolds  number, 

.  .  .  +  + 

as  can  be  inferred  from  the  definition  of  u  ,  u  = 

civ 

The  preceding  comments  are  largely  conjectural 
because  of  the  limited  data,  however,  they  might  provide  an 
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empirical  way  of  examining  and  comparing  the  development  of 
velocity  profiles  in  quantitative  terms,  and  consistent  with 
the  concept  of  fully  developed  profiles. 

The  original  center  line  velocities  measured  with  the 
carriage,  using  the  carriage  static  pressure  probe,  did  not 
yield  the  correct  absolute  value  of  the  velocity  as  stated 
earlier  in  the  "Experimental  Procedure"  section.  However,  the 
deviations  between  these  values  and  the  more  correct  values 
were  found  to  be  constant  for  one  particular  flow  rate,  thus 
making  the  center-line  velocity  dimensionless  by  dividing  the 
value  measured  at  any  position  along  the  tube  axis  by  the  one 
measured  in  the  fully  developed  region,  a  representation  of 
the  center-line  development  was  obtained.  The  center-line 
velocity,  shown  in  Figure  28,  was  measured  for  8  different 
flow  rates,  from  Re  =  8,800  to  55,000.  All  of  the  axial  velo¬ 
city  profiles  which  were  obtained  are  given  in  Tables  38  to  43, 
while  the  data  for  the  two  lower  flow  rates  were  excluded 
from  Figure  28.  At  low  flow  rates,  the  measured  change  in 
impact  pressure  was  small  as  compared  to  the  accuracy  of  the 
micro-manometer  and  the  results  so  obtained  showed  great 
scatter.  The  data  shown  in  Figure  28  covers  a  small  range  of 
Reynolds  numbers,  18,300  to  55,000,  and  the  effect  of  Reynolds 
number  could  only  be  distinguished  for  averaged  center-line 
velocity  values  of  two  and  two  runs  together.  Inspite  of  a 
very  high  center-line  velocity  at  the  entrance,  as  revealed 
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by  Figure  28,  a  consequence  of  jetting  of  the  fluid  from  the 
one-inch  line,  the  velocity  profile  after  a  short  distance 
downstream  had  assumed  the  more  normal  shape  shown  in  Figure 
25.  However,  this  jet  entering  the  6.255  inches  tube,  from  a 
one-inch  pipe  with  a  high  turbulence  leve,  enables  the  flow  to 
attain  a  more  developed  velocity  profile  in  a  shorter  distance 
downstream  than,  for  example,  a  rounded  entrance  such  as  that 
used  by  Deissler(3).  This  would  explain  the  relatively  short 
entrance  length  encountered  in  the  present  investigation. 
Defining  the  entrance  length,  L/D,  as  the  distance  required 
for  the  center-line  velocity  to  reach  within  1.0%  of  the 
constant,  fully  developed  value,  it  was  found  that 

L 

-  ^  40  ,  9000  ^  Re  $  55, 000  (71) 

D 

From  Figure  28  it  can  also  be  seen  that  the  entrance  length 
increases  with  increasing  Reynolds  number.  However,  no 
functional  relationship  could  be  derived  because  of  the  re¬ 
stricted  data  available. 
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V.  CONCLUSIONS 

The  change  in  center-line  velocity  along  the  tube 
axis  was  used  as  the  basis  for  measuring  the  velocity  profile 
development.  The  experimental  data  obtained  with  the  travelling 
carriage  indicated  that  fully  developed  flow  was  reached  at  40 
diameters  downstream  from  the  entrance.  A  slight  increase  in 
entrance  length  with  increasing  Reynolds  number  was  also  ex¬ 
perienced,  However,  no  correct  absolute  value  of  the  center- 
line  velocity  along  the  tube  axis  was  obtained.  This  uncertainty 
partially  resulted  from  the  presence  of  the  carriage,  and  par¬ 
tially  from  errors  associated  with  measurements  of  the  static 
pressure . 

For  future  work  with  the  carriage,  it  would  be  ad¬ 
vantageous  to  extend  the  probes  additionally  beyond  the  car¬ 
riage,  and  at  the  same  time  to  change  the  present  carriage 
static  pressure  probe  to  a  disc  type(24)  static  pressure  probe. 
The  aluminum  6-inch  tube  was  not  always  perfectly  circular  and 
therefore  it  would  improve  the  positioning  of  the  carriage  if 
its  teflon  bearings  could  be  mounted  to  the  carriage-arms  on 
small  springs.  This  would  secure  rigid  positioning  of  the 
carriage  at  all  locations  in  the  tube.  In  further  work,  the 
carriage  could  also  be  used  to  obtain  point  velocities  between 
the  center  line  and  the  wall,  and  in  this  way  give  more  infor¬ 
mation  about  the  velocity  profiles  in  the  developing  section. 

This  is  believed  to  be  particularly  important  near  the  inlet. 
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A  set  of  8  velocity  profiles  under  fully  developed 


conditions  were  measured  in  the  6.255  in.  I.D.  aluminum  test 
tube  at  41  diameters  from  the  entrance.  Various  velocity  pro¬ 
file  relationships  were  tested  and  found  to  correlate  the 
experimental  data  more  or  less  accurately.  Over  the  range  of 
Reynolds  numbers  tested,  Re  =  9,000  to  56,000,  an  equation  based 
on  a  derivation  of  Pai  was  found  to  correlate  the  obtained 
velocity  profiles  very  wells 


u 


(42) 


This  equation  is  more  complicated  than  the  logarithmic  type 
of  distribution,  but  it  is  of  an  easy  form  assuming  a  computer 
is  available.  The  equation  has  the  great  advantage  of  being 
valid  from  the  wall  to  the  center  line.  The  improved  universal 
velocity  distribution  correlation  suggested  by  Bogue,  and 
simplified  slightly  in  the  present  investigation,  was  comparable 
to  the  Pai  correlation,  but  has  the  disadvantage  of  not  being 
valid  near  the  wall.  The  form  of  the  improved  universal 
velocity  distribution  which  was  obtained  in  the  range  of  Rey¬ 
nolds  numbers  from  9,000  to  56,000  was 
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C  (— )  =  1.07  exp  ( 

R 


-(-  -  0 . 77 ) 2 

R 


- - _)  (71) 


0.16 


82 


However,  this  velocity  distribution  could  be  joined  in  a  con¬ 
tinuous  manner  with  one  of  the  equations  derived  for  the  region 
near  the  wall,  if  such  an  extension  were  desirable. 

The  data  obtained  close  to  the  wall  were  not  satis¬ 
factory,  and  the  only  evident  conclusion  was  that  for  y+  <  25 
viscous  effects  become  significant  making  the  universal  velocity 
distribution  no  longer  valid.  This  is  in  agreement  with  Deissler 

who  found  the  universal  velocity  distribution  to  be  valid  for 
+ 

y  ^>26.  The  different  equations  suggested  for  the  wall  region, 
y+  <  26,  still  remain  to  be  tested  experimentally.  However, 
the  present  experimental  equipment  does  not  provide  the  accur¬ 
acy  needed  in  measuring  point  velocity  close  to  the  wall. 

The  pressure  drop  data  for  friction  factor  calcula¬ 
tions  were  obtained  in  the  upper  region  of  the  tube.  These 
results  gave  slightly  higher  friction  factors  than  the  Blasius 
friction  factor  correlation  for  smooth  pipes  would  predict. 
However,  the  region  where  the  pressure  drop  data  were  obtained 
was  found  to  deviate  locally  from  a  circular  cross-section, 
and  thus  the  obtained  friction  factors  were  questionable. 

Further  investigation  of  pressure  drop  along  other  sections  of 
the  tube  might  clarify  the  uncertainty  regarding  the  true  value 
of  the  friction  factor.  A  survey  of  pressure  drop  along  the 
tube  is  restricted  by  the  length  of  the  tube  and  furthermore 
by  the  developing  shear  stress  in  the  region  near  the  inlet. 

At  low  rates,  the  measurement  of  pressure  drop  along  the  tube 
becomes  greatly  restricted  by  the  decreasing  accuracy  of  the 


micromanometer . 
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NOMENCLATURE 

al'  a2 

Variables  in  the  Pai  correlation 

A,  B 

Constants  in  the  Universal  Velocity  Distribution 

c ,  d 

Variables  in  the  Gill  and  Scher  correlation 

C,  C  • 

Constants  of  integration 

C(y, f ) , 

C(y)  Bogue  correction  functions 

D 

Diameter  of  the  tube,  in. 

F 

.  2 

Cross  section  of  the  tube,  m 

f 

Blasius  friction  factor 

9c 

Dimensional  conversion  factor,  lbmf t/lb^sec2 

Ah 

Pressure  difference,  in. 

K 

Universal  constant 

1 

Mixing  length 

L 

Entrance  length,  in. 

M 

Variable  in  the  Pai  correlation 

n 

Power  law  exponent 

n 

Constant  of  proportionality 

P 

2 

Absolute  pressure,  lb/ft 

po 

Absolute  pressure  at  standard 
conditions  (760  mm  Hg  and  273°K) 

AP 

2 

Dynamic  or  impact  pressure,  lb/ft 

Q 

Flow  rate,  ft^/sec. 

Qo 

Flow  rate  at  standard  conditions,  SCFM 

r 

R 


Re 


Radius,  in. 

Radius  of  tube,  in. 
Reynolds  number 


It  ! 
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S 

T 

T 

u 


u 


av 


um 

UFD 

u+ 

u* 


U4;  U5 
X 


Y 


Y 

€ 


M- 

P 

T 

tl 

T  ip 

TW 

0 


Variable  in  the  Pai  correlation 
Temperature,  °R  or  °K 
Standard  temperature,  273°K 
Point  velocity,  ft/sec. 

Average  or  bulk  velocity,  ft/sec. 

Maximum  velocity,  ft/sec. 

Fully  developed  center-line  velocity,  ft/sec. 
Dimensionless  velocity 
Friction  velocity,  ft/sec 

Universal  constants  in  the  Wasan  correlation 
Distance  downstream  from  entrance  of  tube,  in. 
Distance  from  the  wall,  in. 

Dimensionless  distance  from  the  wall 
Expansion  factor 

2 

Coefficient  of  eddy  diffusivity,  ft  /sec. 
Dynamic  viscosity,  lb/ft. sec. 

Kinematic  viscosity,  ft/sec. 

Density,  lb/ft^ 

2 

Shear  stress,  lb/ft 

Shear  stress  produced  by  viscosity,  lb/ft^ 
Shear  stress  produced  by  turbulence,  lb/ft 
Shear  stress  at  the  wall,  lb/ft 
Variable  in  the  Gill  and  Scher  correlation 


’ 
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APPENDIX  A 
EXPERIMENTAL  AND  CALCULATED  DATA 


38 


1 •  Point  Velocity  Calculations  from  Pitot  Tube  Measurements 


The  standard  equation  to  calculate  velocity  from 
measured  impact  pressure  is: 


u  = 


2g  AP 


air 


.-/here 


u 

AP 


air  velocity,  ft/sec 
32.17,  lbmft/lbfsec2 
impact  pressure,  Ib/ft 


Pair  -  air  density,  lb/ft 

The  air  density  was  obtained  by  assuming  the  ideal  gas  law 


holds 


air 


air 


T  P 
_o 

OTP, 


where 


and 


air, 


=  0.0808  lb/ft J  * 


T  =  27  3°K 


P0  =  760  mm  Hg 


o 

T 


=  2  9  6  °  K 


-4, 


air 


=  0.98387  (10  )  P 


*  Perry,  J.H.,  "Chemical  Engineers'  Handbook",  McGraw-Hill 
Book  Company  Inc.,  New  York,  N.Y.,  Fourth  Edition,  1963. 
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The  impact  pressure  was  measured  in  inches  of  alcohol 


where 


AP 


palc 


Ah 

12 


Ah 

palc 

Density  of  the 


=  length  of  column,  in.  ale, 
=  density  of  alcohol,  lb/ft3 

n-Butyl  alcohol* 


pa^c  =  0.8176  -  0.0004  (t°F  -  50),  gr/cm3 

At  7  5°F 

palc  =  0.8076  gr/cm3  =  50.4185  lb/ft3 

At  a  flowing  air  temperature  of  22 °C  and  a  micromanometer 
temperature  of  75°F 

u  =  l\J  2.75  (106)  — 

where 

Ah  =  length  of  column,  in.  ale. 

P  =  pressure  at  probe  location,  mm  Hg 
The  above  relation  was  used  for  all  velocity  calculations. 


* 


Flow  Corporation,  Model  MM3  Micromanometer  Notes  and 
Instructions,  Bulletin  No.  19,  March  1961. 
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2 .  Flow  Rate  Calculations  from  Orifice  Pressure  Drop  Measurements 


The  standard  equation  to  calculate  flow  rate  is 


Q  =  CFY 


2gc  ap 


pair 


Substituting  the  internal  cross  section  of  the  orifice, 

2 


F  = 


7T  6.2  55' 
(4) (144) 


ft' 


and  the  air  density 


p  =  0.98387  (10“4)  P 

air 


the  standard  equation  becomes 


Q  =  0.7612  CY 


APT 


Expressed  in  terms  of  standard  cubic  feet  per  minute, 


where 


Qo 

c 

Y 

AP 

P 

T 


=  31.73  CY 


APP 


T 


flow  rate,  SCFM 

flow  coefficient,  dimensionless 
expansion  factor,  dimensionless 
pressure  difference  across  orifice,  mm  Hg 
pressure  on  upstream  side  of  orifice,  mm  Hg 
temperature,  °R 
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3 .  Friction  Factor  Calculations  from  Pressure  Dr  op  Measurements 


The  friction  factor  for  incompressible  flow  is 
obtained  from  the  Blasius  equation, 


f  = 


2g, 


D 


’air  uav 


dP 

(-  — ) 

dx 


where 


and 


,  dpv  Pale 

(_  — )  ^  — —  Ah 


50.4 


dx 


D 


Ax 


6.255  in . 


(6) (12) 


Ah 


Upon  substitution,  the  friction  factor  equation  becomes 


where 


28.12 


Ah 


D  •  U 

Kair  av 


f  = 
Ah  = 

pair 


Blasius  friction  factor,  dimensionless 
length  of  alcohol  column,  in. 

O 

air  density,  lb/ft 
average  velocity,  ft/sec 
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TABLE  1 

Entrance  Velocity  Profiles  -  Old  Vanes 


Point  of  measurement 
Barometric  pressure 
Flowing  temperature 
Pressure  at  probe  location 
Orifice  pressure 
Orifice  pressure  difference 
Reynolds  number 


4.2  diameter 
708.0  mm  Hg 
531.6  °R 
709.3  mm  Hg 
723.1  mm  Hg 
9 . 64  mm  Hg 
18 , 000 


Probe  Position 

(in.) 


Ah 

(in.  ale.) 


u 

( ft/sec) 


0.00 

0.0123 

6.92 

1.00 

in 

0.0132 

7.17 

2  .00 

in 

0.0125 

6.97 

3.00 

in 

0.0094 

6.05 

1.00 

out 

0.0111 

6 . 56 

2.00 

out 

0.0094 

6.05 

3.00 

out 

0.0061 

4.87 

TABLE  2 


Entrance 

Velocity  Profiles  -  New  Vanes 

Point  of 

r  measurement 

4.2  diameter 

Barometric  pressure 

693.4 

mm 

Hg 

Flowing 

temperature 

531.6° 

R 

Pressure 

:  at  probe  location 

69  5.5 

mm 

Hg 

Orifice 

pressure 

715.2 

mm 

Hg 

Orifice 

pressure  difference 

15.3  mm  ] 

Hg 

Reynolds 

;  number 

23,  000 

Probe  Position 

Ah 

u 

(in.) 

(in.  ale 

.) 

( ft/sec) 

0.00 

0.0216 

9  .25 

1.00 

in 

0.0194 

8.77 

2  .00 

in 

0.0163 

8.04 

2  ,75 

in 

0.0131 

7  .20 

3.00 

in 

0.0085 

5.81 

1.00 

out 

0.0206 

9.04 

2.00 

out 

0.0173 

8.28 

3.00 

out 

0.0119 

6.86 

94 


TABLE  3 

Calibration  of  Orifice 


Re 

Qo 

SCFM 

Q0/C 

SCFM 

C 

3,874 

35.588 

54.944 

0.6477 

13,  117 

52 .607 

81.067 

0.6489 

18, 176 

72 .898 

113.554 

0.6420 

23,  816 

95.516 

144.288 

0.6620 

28,987 

116.255 

176.394 

0.6591 

32,711 

131.191 

209  .238 

0.6270 

38, 181 

153.128 

243.383 

0.6292 

56 , 045 

224.773 

355.982 

0.6314 

TABLE  4 

Friction  Factor 

Re 

Ah 

(in.  ale.) 

P  -3  uav 

(lbs/ft3)  (ft/sec) 

f 

8,  874 

0.0010 

0.06835 

3.058 

0 .04399 

13, 117 

0.0017 

0.06881 

4.49 

0.03446 

18, 176 

0.0031 

0.06917 

6.19 

0.03289 

23, 816 

0.0045 

0.06883 

8.15 

0.02767 

28,987 

0.0063 

0.06953 

9.82 

0.02642 

32,711 

0.0085 

0.07011 

10.99 

0.02822 

38, 181 

0 . 0109 

0.06972 

12  .90 

0.02641 

53, 696 

0.0197 

0.07015 

18.03 

0.02429 

• 
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TABLE  5 

Traverse  at  Right  Angle 


Pt.  of  measurement 
Baro.  press. 

Flowing  temp . 

Press  at  probe  loc . 
Orifice  pressure 
Orifice  press,  diff. 
Reynolds  number 


41  diameter 
702.6  mm  Hg 
531. 6°R 
703.8  mm  Hg 
718.5  mm  Hg 
9 . 56  mm  Hg 
18, 000 


Main  Position  At  90° 


Probe 

Position 

( in) 

Ah 

(  in . ale  . ) 

u 

(  ft) 

sec 

u(  av) 
ft 

( - ) 

sec 

Ah 

(in.  ale.) 

U 

ft 

(  ) 

sec 

u(  av) 
ft. 

(  ) 

sec 

0.00 

0.0150 

7  .66 

7  .66 

0.0151 

7.68 

7.68 

0.50 

in 

0.0145 

7.53 

7 .52 

0.0147 

7.58 

7 . 59 

1.00 

in 

0.0140 

7  .40 

7.36 

0.0142 

7  .45 

7.38 

1.50 

in 

0.0129 

7 .10 

7  .09 

0.0128 

7.07 

7  .01 

2.00 

in 

0.0110 

6.56 

6.59 

0.0112 

6.61 

6.58 

2.50 

in 

0.0094 

6 . 06 

6.06 

0.0092 

5.99 

5.95 

2.75 

in 

0.0078 

5.52 

5.52 

0.0080 

5.59 

5.54 

3.00 

in 

0.0050 

4.42 

4.51 

0.0035 

3.70 

4.20 

3.05 

in 

0.0031 

3.48 

3.79 

0.0010 

1.98 

2 .52 

0.50 

out 

0.0144 

7.50 

0.0147 

7.59 

1.00 

out 

0.0137 

7.32 

0.0136 

7 . 30 

1.50 

out 

0.0128 

7  .07 

0.0123 

6.94 

2.00 

out 

0.0112 

6.61 

0.0109 

6.54 

2,50 

out 

0.0094 

6.06 

0.0089 

5.91 

2.75 

out 

0.0078 

5.52 

0.0077 

5.49 

3.00 

out 

0.0054 

4.59 

0.0056 

4.69 

3.05 

out 

0.0043 

4.10 

0.0040 

3.96 
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TABLE  6 

Traverse  at  Right  Angle 


Pt.  of  measurement 
Baro.  press. 

Flowing  temp . 

Press,  at  probe  loc . 
Orifice  press. 
Orifice  press,  diff. 
Reynolds  number 


4.2  diameter 
705.8  mm  Hg 
531. 6°R 

708.1  mm  Hg 

754.2  mm  Hg 
16.1  mm  Hg 
24, 000 


Main  Position 

At  90° 

Probe 

Ah 

u 

u(av) 

Ah 

u 

u(  av) 

Position 
( in . ) 

( in . ale  . ) 

ft 

( - ) 

sec 

(— ) 
sec 

(in.  ale.) 

ft 

(• - ) 

sec 

£  t 

sec 

0.00 

0.0230 

9.45 

9.45 

0  o  02  31 

9.47 

9.47 

1.00 

in 

0.0221 

9.27 

9.29 

0.219 

9.22 

9.24 

2.00 

in 

0.0198 

8.77 

8.77 

0.0198 

8.77 

8.77 

2.50 

in 

- 

— 

— 

0.0182 

8.41 

8.44 

3.00 

in 

0.0129 

7  .08 

7.15 

— 

— 

- 

1.00 

out 

0.0223 

9.31 

0.0221 

9.26 

2 .00 

out 

0.0198 

8.77 

0.0198 

8.77 

2.50 

out 

— 

— 

0,0185 

8.47 

3.00 

out 

0.0134 

7  .21 

— 
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TABLE  7 

Fully  Developed  Velocity  Profile 


Pt .  of  Measurement 
Baro .  press 
Flowing  temp . 

Press,  at  Probe  loc . 
Density 
Orifice  press 
Orifice  press  diff. 
Expansion  factor 


41  diameter 
695.5  mm  Hg 
531. 6°R 
695.5  mm  Hg 
0.06835  lbs/ft3 
699 .0  mm  Hg 
2.28  mm  Hg 
1.0 


Probe  Position 

Ah 

u 

u(  av) 

ur 

ft 

ft 

f t . in  . 

u 

(in.) 

(in.  ale.) 

(— ) 

( - ) 

( - — ) 

u„ 

sec 

sec 

sec 

m 

0.00 

0.0037 

3.82 

3.82 

0.0 

1.0000 

0.50 

in 

0.0036 

3.77 

3.75 

1  .875 

0.9817 

1.00 

in 

0.0035 

3.72 

3.70 

3.700 

0.9686 

1.50 

in 

0.0033 

3.61 

3.56 

5.340 

0.9319 

2.00 

in 

0.0028 

3.33 

3.27 

6.540 

0.8560 

2.50 

in 

0.0025 

3.14 

3.08 

7  .700 

0.3063 

2.75 

in 

0.0018 

2  .67 

2  .72 

7  .480 

0.7120 

3.00 

in 

0.0009 

1.89 

2.08 

6.240 

0.5445 

0.50 

out 

0.0035 

3.72 

1.00 

out 

0.0034 

3.67 

1.50 

out 

0.0031 

3.50 

2  .00 

out 

0.0026 

3.21 

2.50 

out 

0.0023 

3.02 

2.75 

out 

0.0019 

2  .74 

3.00 

out 

0.0013 

2  .27 

Average  velocity 

3.085  ft/sec 

Flow  rate 

35.588  SCFM 

Flow  coefficient 

0.6477 

Reynolds  numb^ 

8,874 

Ratio  of  av .  to  max.  vel . 

0.80 

Friction  factor 

0.04399 

> 
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TABLE  8 

Fully  Developed  Velocity  Profile 


Pt .  of  measurement 
Baro.  press. 

Flowing  temp . 

Press,  at  probe  loc . 
Density 

Orifice  press. 
Orifice  press,  diff. 
Expansion  factor 


41  diameters 
699.6  mm  Hg 
531. 6°R 

700.2  mm  Hg 
0.06881  lbs/ft3 

707.3  mm  Hg 
4.93  mm  Hg 
.79 


Probe 

( 

Position 

in . ) 

Ah 

(in.  ale.) 

u 

(— ) 
sec 

u(av) 

(— ) 
sec 

ur 

^  f t . in  . 
sec 

u 

um 

m 

0.00 

0.0082 

5.68 

5.68 

0.0 

1.0000 

0.50 

in 

0.0080 

5.60 

5.60 

2.800 

0.9859 

1.00 

in 

0.0076 

5.47 

5.49 

5.490 

0.9665 

1.50 

in 

0.0068 

5.17 

5.19 

7  .785 

0,9137 

2.00 

in 

0.0059 

4.81 

4.77 

9.540 

0.8398 

2.50 

in 

0.0050 

4.43 

4.43 

11.075 

0.7799 

2.75 

in 

0.0043 

4.11 

4.11 

11.3025 

0.7236 

3.00 

in 

0.0020 

2.80 

3.03 

9.090 

0.5335 

3.05 

in 

0.0008 

1.77 

2  .29 

6.9845 

0.4032 

0.50 

out 

0.0080 

5.60 

1.00 

out 

0.0077 

5.50 

1.50 

out 

0.0069 

5.21 

2  .00 

out 

0.0057 

4.73 

2.50 

out 

0.0050 

4.43 

2.75 

out 

0.0043 

4.11 

3.00 

out 

0.0027 

3.26 

3.05 

out 

0.0020 

2.80 

Average  velocity 

4.49  ft/sec 

Flow  rate 

52.607  SCFM 

Flow  coefficient 

0.6489 

Reynolds  number 

13, 117 

Ratio  of  av.  to  max.  velocity  0.79 
Friction  factor  0.03446 
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TABLE  9 

Fully  Developed  Velocity  Profile 


Pt .  of  measurement 
Baro.  press. 

Flowing  temp . 

Press,  of  probe  loc . 
Density 
Orifice  press 
Orifice  press,  diff. 
Expansion  factor 


41  diameter 
702.6  mm  Hg 
531. 6°R 
7  0  3 . 8  mm  Hg 
0.06917  lbs/ft3 
718.5  mm  Hg 
9 . 56  mm  Hg 
0.9955 


Probe 

Position 

Ah 

u 

u(  av) 

ur 

( 

in . ) 

(in.  ale.) 

(— ) 

sec 

(— > 
sec 

f t . in  .  ^ 

sec 

u 

um 

0.00 

0.0150 

7  .66 

7  .66 

0.0 

1.0000 

0.50 

in 

0.0145 

7.53 

7 .52 

3.760 

0.9817 

1.00 

in 

0.0140 

7  .40 

7.36 

7.360 

0.9608 

1.50 

in 

0.0129 

7.10 

7  .09 

10.635 

0.9256 

2.00 

in 

0.0110 

6 . 56 

6.59 

13.180 

0.8603 

2.50 

in 

0.0094 

6 . 06 

6 . 06 

15.150 

0.7911 

2.75 

in 

0.0078 

5.52 

5.52 

15.180 

0.7206 

3.00 

in 

0.0050 

4.42 

4.51 

13.530 

0.5888 

3.05 

in 

0.0031 

3.48 

3.79 

11.5595 

0.4948 

0.50 

out 

0.0144 

7.50 

1.00 

out 

0.0137 

7 . 32 

1.50 

out 

0.0128 

7  .07 

2.00 

out 

0.0112 

6.61 

2.50 

out 

0.0094 

6.06 

2.75 

out 

0.0078 

5.52 

3.00 

out 

0.0054 

4.59 

3.05 

out 

0.0043 

4.10 

Average  velocity 

6.19  ft/sec 

Flow  rate 

72.898  SCFM 

Flow  coefficient 

0.6420 

Reynolds  number 

18, 176 

Ratio  of  Av.  to  max. 

velocity  0.808 

Friction  factor 

0.03289 
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TABLE  10 

Fully  Developed  Vel ocity  Pro file 


Pt .  of  measurement 
Baro .  press. 

Flowing  temp. 

Press,  at  probe  loc . 
Density 

Orifice  press . 
Orifice  press,  diff. 
Expansion  factor 


41  diameter 
698.2  mm  Hg 
531. 6°R 
700.4  mm  Hg 
0.06883  lbs/ft3 
722.7  mm  Hg 
15.44  mm  Hg 
0.9925 


Probe 

( 

Position 

in  . ) 

Ah 

(in.  ale.) 

u 

ft 

( — > 

sec 

u(  av) 

(— ) 
sec 

ur 

f t . in  . 

\  7 

sec 

u 

m 

0.00 

0.0256 

10.03 

10.03 

0.000 

1.0000 

0.50 

in 

0.0250 

9.91 

9.94 

4.970 

0.9910 

1.00 

in 

0.0234 

9.59 

9.65 

9.650 

0.9621 

1.50 

in 

0.0213 

9.15 

9.22 

13.830 

0.9192 

2.00 

in 

0.0193 

8.70 

8.69 

17.380 

0.8664 

2 ,50 

in 

0.0161 

7  .95 

7  .99 

19 .975 

0.7966 

2.75 

in 

0.0136 

7.31 

7 .36 

20.240 

0.7338 

3.00 

in 

0.0077 

5.50 

5.92 

17  .760 

0.5902 

3,05 

in 

0.0032 

3.55 

4.74 

14.457 

0.4726 

0.50 

out 

0.0253 

9.97 

1.00 

out 

0.0240 

9  .71 

1.50 

out 

0.0220 

9.29 

2.00 

out 

0.0192 

8.68 

2.50 

out 

0.0164 

8.02 

2.75 

out 

0.0140 

7  .41 

3.00 

out 

0.0102 

6.33 

3.05 

out 

0.0090 

5.93 

Average  velocity 

8.15  ft/sec 

Flow  rate 

95.516  SCFM 

Flow  coefficient 

0.6620 

Reynolds  number 

23, 816 

Ratio  of  average  to  max.  valocity 

0.813 

Friction  factor 

0.02767 
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TABLE  11 

Fully  Developed  Velocity  Profile 


Pt«  of  measurement 
Baro .  press. 

Flowing  temp . 

Press,  at  probe  loc . 
Density 
Orifice  press. 
Orifice  press,  diff. 
Expansion  factor 


41  diameter 
704.3  mm  Hg 
531. 6°R 
7  07 . 5  mm  Hg 
0.06953  lbs/ft3 
740.0  mm  Hg 
22.65  mm  Hg . 
0.99 


Probe  Position 
(in.) 


Ah 

(in.  ale.) 


u 

(— ) 

sec 


u(  av) 
ft 

( - ) 

sec 


ur 


^  f t . in  . 


sec 


u 


u 


rn 


0.00 

0.0334 

12  .22 

12  .22 

0.000 

1.0000 

0.50 

in 

0.0375 

12.07 

12.10 

6.050 

0.9902 

1.00 

in 

0.0356 

11.76 

11.83 

11.830 

0.9681 

1.50 

in 

0.0327 

11.27 

11.35 

17  .025 

0.9288 

2  .00 

in 

0.0283 

10.49 

10.54 

21.080 

0.8625 

2.50 

in 

0.0231 

9  *48 

9.53 

23.825 

0.7799 

2  .75 

in 

0.0199 

8.79 

8.88 

24.420 

0.7267 

3.00 

in 

0.0124 

6.94 

7.11 

21.330 

0.5818 

3.05 

in 

0.0052 

4.50 

5.14 

15.677 

0.4206 

0.50 

out 

0.0378 

12.12 

1.00 

out 

0.0364 

11.89 

1.50 

out 

0.0336 

11.43 

2.00 

out 

0.0288 

10.58 

2.50 

out 

0.0236 

9.58 

2.75 

out 

0.0207 

8.97 

3.00 

out 

0.0139 

7  .27 

3.05 

out 

0.0086 

5.78 

Average  velocity 

9.82  ft/sec 

Flow  rate 

116.255  SCFM 

Flow  coefficient 

0.6591 

Reynolds  number 

28,987 

Ratio  of  average  to  max.  velocity 

0.804 

Friction  factor 

0.02642 
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TABLE  12 

Fully  Developed  Velocity  Profile 


Pt.  of  measurement 
Baro.  press. 

Flowing  temp . 

Press,  at  probe  loc. 
Density 
Orifice  press. 
Orifice  press,  diff. 
Expansion  factor 


41  diameter 
709.0  mm  Hg 
531 .6°R 
713.4  mm  Hg 
0.07011 
760.0  mm  Hg 
31 .2 5  mm  Hg 
0.9865 


Probe  Position  Ah 

(in.)  (in.  ale.) 


0.00 

0.0471 

0.50 

in 

0.0461 

1.00 

in 

0.0440 

1.50 

in 

0.0403 

2.00 

in 

0.0359 

2.50 

in 

0.0296 

2.75 

in 

0.0253 

3.00 

in 

0.0161 

3.05 

in 

0.0100 

0 , 50 

out 

0.0461 

1.00 

out 

0.0439 

1.50 

out 

0.0404 

2.00 

out 

0.0355 

2,50 

out 

0.0301 

2.75 

out 

0.0264 

3.00 

out 

0.0187 

3.05 

out 

0.0125 

Average  velocity 
Flow  rate 
Flow  coefficient 
Reynolds  number 
Ratio  of  average  to 
Friction  factor 


u  u(av)  ur 

ft  ft  ft. in. 

(— )  (~ ~)  (— ) 

sec  sec  sec 


13.47 

13.47 

0.0 

1,0000 

13.33 

13,33 

6.667 

0.9896 

13.02 

13.02 

13.020 

0.9666 

12  .46 

12.47 

18.705 

0  .92  58 

11 .76 

11.73 

23.460 

0.8708 

10.68 

10,73 

26,825 

0.7966 

9.88 

9  .98 

27 .445 

0.7409 

7  .88 

8.19 

24.570 

0,6080 

6.21 

6.53 

20.069 

0.4885 

13.33 
13.01 
12  .48 
11.70 
10,77 
10.03 
8.49 
6,94 


10.99  ft/sec 
131.191  SCFM 
0.6270 
32,711 

max.  velocity  0.02822 

0.816 
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TABLE  13 

Fully  Developed  Velocity  Profile 


Pt.  of  measurement 
Baro.  press. 

Flowing  temp. 

Press,  at  probe  loc . 
Density 
Orifice  press. 
Orifice  press,  diff. 
Expansion  factor 


41  diameter 
703.0  mm  Hg 
531. 6°R 
709.4  ram  Hg 
0,06972  lbs/ft3 
774.0  mm  Hg 
41,9  mm  Hg 
0.932 


Probe 

( 

Position 

in . ) 

Ah 

(in.  ale.) 

u 

ft 

{  \ 

u(av) 

ft 

/  \ 

ur 

f t . in  . 

/  ... _ _ \ 

u 

u 

m 

\  ) 
sec 

\  ) 
sec 

\  ) 
sec 

0.00 

0.0636 

15.70 

15.70 

0.0 

1.0000 

0.50 

in 

0.0625 

15.57 

15.61 

7  .805 

0.9943 

1.00 

in 

0.0599 

15.24 

15.24 

15.240 

0.9707 

1.50 

in 

0.0550 

14.64 

14.62 

21.930 

0,9312 

2,00 

in 

0.0482 

13.67 

13.71 

27.420 

0.87  32 

2.50 

in 

0.0405 

12 .53 

12 .52 

31,300 

0.7975 

2.75 

in 

0.0350 

11.65 

11.63 

31.983 

0.7408 

3.00 

in 

0.0262 

10.03 

10.02 

30.060 

0.6382 

3.05 

in 

0.0222 

9 ,28 

9  .28 

28.304 

0.5911 

0.50 

out 

0.0631 

15.64 

1.00 

out 

0.0599 

15.24 

1.50 

out 

0.0549 

14.59 

2  .00 

out 

0.0488 

13.75 

2.50 

out 

0.0404 

12.51 

2.75 

out 

0.0347 

11.60 

3.00 

out 

0.0256 

9  .96 

3.05 

out 

0.0222 

9.28 

Average  velocity 

12.90  ft/sec 

Flow 

rate 

153,128  SCFM 

Flow 

coefficient 

0.6292 

Reynolds  number 

38, 181 

Ratio  of  average  to  max,  velocity 

0.822 

Friction  factor 

0,02641 
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TABLE  14 

Fully  Developed  Velocity  Profile 


Pt.  of  measurement 
Baro .  press. 

Flowing  temp . 

Press,  at  probe  loc . 
Density 
Orifice  press. 
Orifice  press,  diff. 
Expansion  factor 


41  diameter 
699 .0  mm  Hg 
531. 6°R 
713.0  mm  Hg 
0.07009  lbs/ft3 
952.0  mm  Hg 
74.0  mm  Hg 
0.9745 


Probe  Position 

( in . ) 


Ah 

(in.  ale.) 


u 
ft 


( - ) 

sec 


u(  av) 
ft 

( - ) 

sec 


ur 

ft .  in 

( - _) 

sec 


u 


0.00 

0.1337 

0.50 

in 

0.1311 

1.00 

in 

0.1243 

1.50 

in 

0.1151 

2.00 

in 

0.1034 

2.50 

in 

0.0867 

2.75 

in 

0.0755 

3.00 

in 

0.0563 

3.05 

in 

0.0487 

0.50 

out 

0.1315 

1.00 

out 

0.1262 

1.50 

out 

0.1173 

2  .00 

out 

0.1034 

2.50 

out 

0.0861 

2.75 

out 

0.0744 

3.00 

out 

0.0554 

3.05 

out 

0.0473 

Average  velocity 
Flow  rate 
Flow  coefficient 
Reynolds  number 
Ratio  of  average  to 


22.71 

22  .71 

0.0 

1.0000 

22  .49 

22.51 

11.255 

0.9912 

21.90 

21.98 

21.980 

0.9679 

21.07 

21.17 

31.755 

0.9322 

19  .97 

19  .97 

39 .940 

0.8793 

18.28 

18.25 

45.625 

0.8036 

17.06 

17  .00 

46.750 

0.7486 

14.73 

14.68 

44.040 

0 . 6464 

13.71 

13.61 

41.511 

0.5993 

22 .52 
22.06 
21.27 
19  .97 
18.22 
16.94 
14.62 
13.51 


18.84  ft/sec 
224.773  SCFM 
0.6314 
56,045 

max.  velocity  0.830 


r 

00.  I 
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4.  Power  Lav/  Distribution 


The  power  law  distribution  is  of  the  following 


form: 


u 

u 


m 


<*) 1/n 

R 


(1) 


The  experimental  data  were  correlated  by  using  three  sets  of 
n-values.  The  first  two  sets  were  obtained  by  the  method  of 
least  squares,  as  follows: 

9 

y  iV 

L)  R 

i) 


ii) 


n  = 


y  u 

In  (— )  In  ( — ) 
R  u 


m 


8 


(-)2 

R 


n  = 


8 


Y  u 

In  (— )  In  ( — ) 
R  u 

m 


The  third  set  of  n-values  was  obtained  from  equation  (61) 


iii) 


n  = 


1 

7 


Experimental  data  and  calculated  values  are  given  in  Tables 


15  to  18. 
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TABLE  19 


Re 

Power 

Law  Exponent-Velocity 

Ratio 

Equation 

(43) 

uav 

um 

Measured 

uav 

n  — — 

um 

Equation 

(61) 

n 

Equation 

(63) 

uav 

um 

8,  874 

5.863 

0.800 

5.542 

0.777 

0.766 

13, 117 

4.863 

0.790 

5.820 

0.786 

0.779 

18, 176 

5.834 

0.808 

6.062 

0.793 

0.789 

23,  816 

5.724 

0.813 

6.270 

0.799 

0.796 

28,987 

5.239 

0.804 

6,  426 

0.803 

0.801 

32,711 

5.986 

0.816 

6.524 

0.805 

0.804 

38, 181 

7  .124 

0.822 

6.651 

0.809 

0.807 

56, 045 

7 .337 

0.830 

6.978 

0.816 

0.816 

Ill 


5  o  Semilogar ithmic  Velocity  Distribution 


i) 

um  ~  u  1 

In 

(-) 

(8) 

u*  I< 

— 
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um  -  u  1 

In 

(1  -A 

/i  -  -  )  +  A 

1  -1 

(12) 

u*  I< 

\ 

/  R  \| 

R 

The  constant, . K,  was  assumed  to  be  universal  and  was  calculated 
by  the  method  of  least  squares  for  each  of  the  two  equations, 
for  equation  8  n 

V"  y  2 

)  (In  (-)  ) 

Lx  R 


n 


z 


um  ~  u  y 

— -  (In  (-)) 

u*  R 


for  equation  12 


Experimental  data  and  calculated  values  are  given  in  Tables 


20  to  23. 
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6  .  Universal  Velocity  Distribution 
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(64) 


The  constants  A  and  B  were  assumed  to  be  universal  and  were 


obtained  by  the  method  of  least  squares 
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(70) 


The  constants  A'  and  B'  were  obtained  as  above,  while  the 
correction  function,  C(y/R) ,  was  determined  by  a  simple  trial- 
and-error  procedure.  The  criterion  used  for  a  best  fit  was 
the  minimum  of  the  sum  of  the  absolute  values  of  percent  error. 
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7 .  Velocity  Distribution  Close  to  the  Wall 


Derivation  of  equation  (65) 


a) 


u+  =  A  In  y+  +  B 


+ 


equation  (64),  valid  at 


at  distance  from  the  wall. 

.  -f  -f  4-4  _L 

b)  u  =  y  4-  Cy  +  Dy 


equation  (32),  valid 


close  to  the  wall 


4"  + 

A  smooth  transition  should  occur  at  y  =  y^ 
Equation  a  gives  for  the  first  and  second  derivatives 


c) 


du' 


dy* 


A 


+ 


d) 


d2u+ 


dv+ 


A 


y 


+‘ 


Equation  b  gives  for  the  first  and  second  derivatives 


e) 


du‘ 

dy 


+ 


-i- 3 

=  1  +  4C  y  +  5D  y 


+ 


f ) 


d2u+ 


=  12C  y  +  20D  y* 


.3 


dy 


.+ 


Equate  c  and  e  and  solve  for  C, 


A 


g) 


c  = 


4y+ 
'  1 


-  -  D  Yl 


4 


y. 


4.3 


Equate  d  and  f  and  solve  for  C, 


A 


h) 


20  + 

c  =  - - -  -  -  D  yi 

12y|  12 


12  3 


Equate  g  and  h  and  solve  for  D, 


i) 


D  = 


A 


Y 


.+ 


5  y 


+  5 


Substituting  expression  for  D  back  into  equation  h  and  obtain 
for  C, 


j) 


C 


15 

12 


Finally,  equating  equation  a  and  b  and  substituting  for  C  and 
D  from  equations  i  and  j,  an  equation  with  only  yj^  as  the  un¬ 
known  is  obtained: 

3  1  B  9 

k)  In  y\  -  -  -  y,  - - +  — 

1  5  A  1  A  20 


Based  on  the  presently  obtained  constants  A  and  B, 

1)  In  y*  -  0.217  y*  =  -  0.994 

By  trial-and-error  solution, 

y+  ^  18 
1 

the  point  at  which  a  smooth  transition  should  occur.  The  value 

of  the  constants  are 

C  =  -  1.403  (10-4) 

D  =  4.55  (10-6) 

The  correlation  valid  for  y'  <  18  is 

u+  =  y4  -  1.403  (10'4)(y+)4  +  4.55  (10"6)(y+)5 

Calculated  values  are  given  in  Table  28. 
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TABLE  28 

Velocity  Distribution  Close  to  the  Wall 


Equation  Equation  Equation 

(  31)  (  33)  (  35) 


1 

1 

y 

+ 

u 

+ 

y 

u+ 

+ 
u  ' 

1.8607 

1 .835 

1.0 

1.000 

1.000 

2 .8360 

2  .753 

2.0 

1 .998 

1 .998 

4.9897 

4.588 

4.0 

3.976 

3.969 

9.1236 

7 . 340 

5.0 

4.944 

4.927 

12 .9388 

9.17  5 

7.0 

6.801 

6.740 

18.1833 

11.010 

10.0 

9.260 

9.052 

21.5816 

11.930 

15.0 

12 .036 

11.352 

25.6905 

12 .840 

20.0 

13.056 

12.112 

30.6900 

13.760 

25.0 

13.965 

14.629 

33.5933 

14.220 
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8 .  Universal  Velocity  Distribution  of  Gill  and  Scher 

The  numerical  integration  of  equation  (39)  was  per¬ 
formed  with  a  4^  order  Legendre-Gauss  quadrature  formula, 

th 

which  is  exact  for  a  7U  degree  polynomial.  The  computer  pro¬ 
gram  used  to  obtain  the  values  given  in  Tables  29  to  32  is 
shown  in  the  next  two  pages.  Changing  the  limit  of  convergence 
from  1%  to  0.1%  barely  changed  the  computed  velocity  profile. 
The  time  required  in  the  latter  case  was  6  minutes  and  12  sec¬ 
onds.  A  numerical  integration  using  the  trapezoidal  rule  was 
initially  attempted  but  it  required  more  than  50,000  increments 
and  more  than  1/2  hour  of  computing  time  to  obtain  the  same 
convergence,  hence  it  was  discarded. 


C 0  =  ( F*XK* ( l.-EXP t-PF*XK/YM) ) ) **2 
H= 1 ,-XK/YM 

TF  ( D.GT.0.0 )  GO  TO  131 
XI =0.0 
GO  TO  1 32 

131  XI  =  ( -1  . +SQRT ( 1 . +4. *GC*n )  ) / 2 • / C  C 

132  X<=H/2,*(  1.+.R61  136311 5 )+A 
CC=(F*XK*( l.-EXP(-PF*XK/YM) ) )  **2 
0= 1 ,-XK/YM 

IF  ( D • G i • 0 • 0 )  GO  i 0  133 
X2  =  0 • 0 
GO  TO  134 

133  X2= ( -1 .+SQRT ( 1 .+4.*CC*D) )/2./CG 

134  XK=H/2.*( 1. -.3399810435 )+A 

CC= ( F*XK* ( 1 ,-FXP (-PF*XK/YM ) ) ) **2 
0= 1 ,-xk/ym 

IF  (n.G'r.0.0)  GO  TO  135 
X  3  =  o .  0 
GO  TO  136 

1  3  5  X3  = ( -1 • +SQRT ( 1  ,+4.*CC*i) )  )/2./CC 

136  XK=H/2.*( l.  +  . 3399810435 ) +A 

CG= (F*XK*( l.-EXP (-PF*XK/YM) ) )**2 
0= 1 ,-XK/YM 

IF  (D.GT.O.O)  GO  TO  137 

X4=0.0 

GO  TO  138 

137  X4= ( -1 • +SORT ( 1 • +4. *CC*D ) ) /2./CC 

138  SI'M=H/?#*(.3478548451*(XI+X2)+.6521451548*(X3+X4) ) 
UC=Uf4-$UM 

A  =  A  +  H 

IF  ( K . L T . N )  GO  TO  13  0 

nFL  =  APS( ur-upr (  i ,  j  )  )  *i 00. /UC 

I F  ( DFL.LT. . 1 )  GO  To  103 

UPC  (  I  ,  J  )  =  UC 

N=  N*2 

FM=FN*2 • 

GO  TO  101 
1  03  UPC (  T ,  J ) =UC 

F  =  UPC(  T  » J ) —UP (  I »  J ) 

FP  =  1 00 • *E /UP (  T  .  J  ) 

FPPOP=FPPOP+ARS(FP) 

WRTTp  (A, 7)  J*N*H*YP(  T  *  J  )  *  UP  (  I  » J  )  *  UPC  (  I  *  J  )  *E*FP»FrRROR 
T  F  (  T . PQ. 8 )  GO  TO  109 
IF  ( J • L  T . 9 )  GO  TO  104 
GO  TO  1 05 

109  IF  (J.LT.8)  GO  TO  104 
CALL  EXIT 
FNfl 


m (\ yx-0 r  *n 

r  r  rs  n  ■  <  ,  ,  *>  #  -  J  3  T 

.  *  r  x 

<  1  r  r»T 

.  '  :  <  » '  .  *>  +  .  r  )  -*■  •  f  —  )  ~  f  x 

M  ■  *  \MrV< 

v  >>  .  )  ti  K  -  .  f  )  *NX* 

M) \  •*-. r  =n 

(  w  «  •  •  O  ) 

.<  =  <v 

V  .  \  (  (  ’  ,  w,(|  '  +  •  f  -  1  =  V 

*>  I  M  f  )  *  .  S  \Ha  NX 

v«  (((  V  »  I  -  !  '  V  -  .  }  X  *  1  )  =  ">*> 

»  y  \  y  .  r  =  i 

r  r  T  o.  (  .  .  -''«<■»)  -u 

.  =  ix 


V  •  \  (  (  •*  *  •  J+  • 

■M  “«•(.)!  •U  ' 
MU  \  >  X  -  -H< 


\.S\  (  ('•*  *  .A  + 

-  •'.  +  (•  •  '  '  i  *  1  A  °  • 


■ 


.  .  (  .  )  -  .  M  .  T  )  C 


ur  -  o:  +  .  [  — )  *  f  x 

+  .  I  )  *  .<s\H=>X 

.  n  ix*  )  =  ■>*> 
m  r  \  >  >  - .  - 
•  v.-  •  .ft)  -)  I 

.  =  kX 

RP  r  it  (  o 

r  )  vi,  + .  r-)r^K 

.  )  *  .  9  \Hs’  I, 

' 11  t  *  +  "M  I  -  )U 

n+-'  =' 

(  / .  r  i .  v )  i  t 

.  f  1  ♦  '  •’  )  —  )  -  M'  1 

r 

')'  «*{  l  .1)7^  ) 
S*  f  =  /i 
.  V  *^s'1  • 

^  '=  (  i  ,  T  )  -•  "H  I 
n  *  T  )  q!  i-  M  ,  T  )  ->Q»  I  =  “> 

. 

( q  )P  4-  ’on  _or»oc- 


. T  >  ■  <  < 

.  *  » 

•  %  ) 

T  c 

or 

rv 

(  « 

•  ^  1 

.  T  ) 

(  ° 

.  • 

) 

.T  i 

•  U 

T  I  X  t  1  I  A  "> 


GILL  AND  SCHFR 


LFGFNDRF  GA'JEF  INTFGRATION 


DIMENSION  V  I ( 9 )  , p ( 9 )  , UP ( 9  *  9 )  >  U ( 9  »  9 ) *  UPC ( 9  »  9 )  »  US ( 9 )  *  UA ( 9 )  *  Y ( 9  *  9 )  *  YP 
1 ( 9  *9)  >RF( 9  )  *BF(9 ) 

1  FORMAT  (1X*9F7.4) 

2  FORMAT  ( 1 X*PF7. 1 ) 


3  FORMAT  (1X*8F6.3) 

4  FORMAT  ( 1X*8F6.1 ) 

5  FORMAT  ( 1HL1HI *  3X , 5HRF ( I  )  > 12X  *2HYM  >9X  * 1HJ  »6X* ] HN*7X*1HH*7X  *7HYP ( I , 
I  J  )  *  3X  ♦  7HUP  (T*J)  >  2X  *  8HUPG  (TO)) 

A  FORMAT  ( 1HL1X* 12 *2X*FB. 1 *?X*F1  0.4 ) 

7  FORMAT  (1H  34X *  12 *1 X  * !6* t X ,F8. 4,2X*F9. ?,?X*F7.3 ,rX*FP ,4*F9. F .Fo.n , 

1  FI  0.5 ) 


READ  (5*1) 
READ  (5*1) 
READ  (5*2) 
READ  (5*3) 
R  F  A  D  (5*4) 


(  (  U  (  10)0=1*9)  *1  =  1*8) 
(!Y(I,J) » J  =  1 *9 )  ,  1  =  1  *8  ) 
( RE (  I  )  *  1  =  1  ,8  ) 

( UA(  I  )  *  1  =  1  *8  ) 

(P(  I  )  *  1  =  1  *8) 


C  =  .  3  1 6 
F  =  .  36 


F  F  =  2  ^  * 

WR  T  TF  (6*5) 


1=0 


31  1=1+1 

RF  (  I  )  =C/RF (  I  ) **. 25 

US  (  I  )  =  UA  (  I  )  OORT  (  PF  (  I  )  / 8  •  ) 

VI  (  I  )  =  • 124933/P (  I  ) 

J=  0 

30  JrJ+1 

UP (  I  * J ) =U(  I  » J ) /US (  I  ) 

YM  I  *  J  )  =Y  (  I  ,  J  )*US  (I)/VI(I)/12. 
IF  (I .FQ.8 )  GO  TO  32 
TF  (J.LT.9)  GO  TO  30 
GO  TO  31 

32  IF  (J.LT.8)  GO  TO  30 
140  1=0 

105  1=1+1 

YM  =  RF (  I  )  / 2  •  *SQR  T ( RF (  I  )/8. ) 
WRITE  ( 6*6 )  I  *  RE (  I  )  * YM 


J=  0 

FRROR=0.0 
104  J=J+1 


N=  10 
FN=10. 

UPC (  I  * J  )  =0#  0 
1  Ol  H  =  YP (  I  * J ) /EN 
UC=0.0 

PF=(  YM  — ^,0  .  )  /  FF 


SUM=rO.  0 
A  =  0  •  0 
K  =  0 

130  K  =  K  + 1 

XK=H/2.*( 1.-.8611363115 ) + A 


'  r  n  ’  *  i  j  t 


.  !  .  (  '  ♦  { 

o  •  (  o  • 

1  !/  1 

•  ( 

«<•>*>}  f 

(  p 

• 

r  ip* 

f 

( 

• 

%■ 

<r 

► 

r ) 

( f 

• 

♦ 

T  A  K’  rf  , 

f 

• 

. 

X 

r ) 

TAvtjn  i 

* 

. 7 ) 

,  '  r .  •  •  .  "  .  1  '  r  . 

,  .  •  .  r  ,  (  I  -  ■<  . 

i  H  r 

■ ) 

I  M..  •  T  '  *  .  •  c  .  ( •  ,  r  )  ^  •  •  ,  *  [  l_  r 


i  .  •  *  *  .  .  » c  .  mm  r  i  r  r  v 

.  .  '  .  .  •  -  .  ^  .  .  ■  .Of 

nr)  T  m  fl  0 r' 

r 

(  ■>.  r*1 

r 

r  .  -J.O*  -  L  »  (  L «  I  )  U 

* 

{  .  -■  T  .  (  f  - 1  ♦  ( l  f  )  ’ )  ) 

(  •  ' 

<  ' .  r  = ! .  (  i  )  ! 

( <;♦  1 )  n.\  sn 

(  R  •  [  -  ♦  (  i  )  ) 

. 

(  Q  t  r  =  T  ,  (  ,  )  4 ) 

< *  ) 

:  f .  * : 

*  .  =0 

.  -- 

!  *  a  )  1  r 

r  +  I  =  1 

.  *  r 

)  >o  -  (  t  i  r 

i .  -  \  (  i  )  )  no.  - ( 1  ) '  -(  !  j  -v 

(  i  )  n  \  ; 

4 s  r . * (  1  1  tv 

l  =L 

f  +  !  si 

(  1  )  \ ( L* 

i  )  -  C  L  *  i  )  n 

•  &f  V  1  )  :  • 

.  )  '(IMH 

<°.0^.T ) 

<  .  1 .  U 

♦ 

r  nT 

(  ft  •  ' .  . 

-  T 

'  A 

(  .  \  (  !  Ma  )  Tno  * 

.  C  \  (  )-jgSMT 

w  *  (  I i  n*  i 

. 

.  rt'r.pq 

4-t  =L  AJ 

o  r 

.  r  =  n  ■» 

.  = ( L* T  ) 

1  ■  o  \  n  ,  i  )  q  >  =  r 

-  •  ’ 

•  — 

\  (  .  a  —m  Y  )  =  q 

.  =V| '? 

C.  =  A 

=  Y 

-  (  A-  )  *  .  ^\H-  v 


Velocity  Distribution  of  Gill  and  Scher 


128 


£ 

oo  (j\ 

CM 

CO 

04 

in 

in 

CO 

CO 

uo 

o 

CO 

oo 

m 

O' 

cm  y 

0 

00  04 

CM 

in  y 

VD 

CO 

oo 

m 

in 

in 

00 

03 

03 

c — 1 

00  O 

•H 

CM  CM 

coy  CNrH 

i — 1 

CO 

ID  CM 

oo 

p" 

CM 

00 

r" 

o  00 

-P 

00  O 

O  03 

CO 

03 

00 

O' 

y  cm 

y 

OM 

03 

vD 

CO 

OO  UO 

£ 

04  O' 

O  03 

00 

vD 

in 

vD 

o  in 

rH 

y 

o 

CM 

oo 

co  in 

> 

f\J  04 

rH 

00 

o 

t — ! 

i — 1 

CM 

o 

o 

rH 

CM 

00 

1 — 1 

y 

y  uo 

0 

1 - 1 

1 

i 

1 

1 

I 

1 

i 

i  i 

P 

1 

£ 

CT>  m 

00 

CO 

m 

y 

i — I 

03 

03 

r-' 

O 

00 

y 

CO 

rH 

VD  CM 

0 

O'  O'* 

O' 

y 

O' 

in 

CO 

vD 

CO 

03 

O 

y 

p" 

y 

CO 

[■"  VD 

•H 

m  i — i 

y 

i — 1 

O' 

y 

ON 

in 

03 

y 

y 

00 

CM 

CO 

uo 

OO  vD 

40 

O' 

+ 

00  oo 

cr> 

y 

i — 1 

vD 

03 

o 

00 

o 

CM 

vD 

O 

r" 

CM 

m  vD 

(0 

oo 

£ 

£ 

03  y 

in 

o- 

CO 

CO 

CO 

03 

00 

rH 

in 

vD 

CO 

CO 

03 

03  03 

i — l 

i — 1 

» — l 

1 — 1 

rH 

rH 

I — I 

rH 

i — 1 

1 - 1 

rH 

rH 

! — 1 

rH  i — 1 

H 

TO 

0) 

04  00 

CO 

CM 

03 

VD 

oo 

rH 

LO 

CO 

t" 

O 

vO 

vD 

vD 

03  CM 

£ 

vd  y 

CO 

vD 

y 

vD 

CM 

00 

03 

o 

vD 

y 

CO 

CM 

CM 

CM  OM 

£ 

+ 

VD  03 

O' 

O* 

04 

03 

OM 

m 

00 

rH 

O 

CM 

y 

o 

1 - 1 

m  co 

W 

£ 

fO 

O  00 

in  *x> 

CO 

CO 

03 

03 

CO  rH 

in  vd 

0"  03 

O 

o  o 

(l) 

i — 1  i — 1 

» — 1 

rH 

rH 

1 - \ 

rH 

rH 

rH 

rH 

rH 

» — 1 

rH 

CM 

CM  CM 

g 

oo  m 

O  i — ! 

00 

y 

in 

VD 

y  y 

in 

in 

t'' 

CD 

o 

rH  00 

00  vO 

cr>  y 

03 

y 

03  y 

t"' 

y 

03  y  y  y  lo 

uo  uo 

+ 

LO  r-4 

o* 

O  04 

in 

O' 

O 

CO 

CM 

O  03 

VD 

OO 

o 

r"  y 

>i 

i — i  y 

VO 

PM 

O' 

CM 

O' 

00 

03 

ID 

CO  03 

00 

O'  rH  y  oo 

H  00 

in 

o  y 

03 

oo  CO 

rH  y 

r"  y 

o 

r*' 

00  03 

rH 

1 — 1 

« — 1 

CM 

CM 

rH 

CM 

CM 

OO  00 

in 

■P 

£ 

m  <d 
o  £ 
<u 

•  p 

o  u 
£  £ 
H 


!Z! 


oooooooo 

VDOOCMPMCMyyy 
rH  CM 


OOOOOOOOO 
i — I  vO 


4J 

£  • 

•H  0 
0  S5 

n. 


rHCMOoymvDr'OO  rHcMooyiniDr-'COcri 


O 

"H  *'“** 

•P  VO  +  S 

03  VO  >1 

£  — 

W 


y 

vO 

in 

o 

oo 

CO 

04 


03 

03 

y 

y 

• 

CO 

03 

oo 


y 

O' 

cc 

00 


O' 


co 


(D 

05 


' 


Velocity  Distri bution  of  Gill  and  Scher 


129 


G 

0 

0 

1 - 1 

P* 

CN 

ro 

00 

p* 

rH 

0 

1 - 1 

y 

0 

0 

rH 

0 

0 

ro 

y 

in 

CN 

CN 

CN 

00 

'st1 

0 

0 

0 

0 

0 

0 

y 

0 

0 

ro 

ro 

1 - 1 

CP 

00 

CN 

0 

0 

y 

O 

0 

0 

P- 

0 

0 

-p 

1 — 1 

y 

m 

0 

CO 

00 

rH 

rH 

0 

0 

P» 

0 

o 

y 

0 

y 

ro 

0 

O 

« — 1 

P> 

O 

<N 

0 

0 

y 

0 

0 

0 

0 

0 

0 

o 

e 

• 

0 

> 

in 

CM 

y 

CN 

CN 

rH 

CN 

0 

y 

0 

.H 

0 

0 

rH 

O 

0 

CD 

1 

1 

l 

1 

i 

1 

1 

i 

Q 

G 

y 

m 

ro 

in 

CN 

0 

1 - 1 

P* 

0 

0 

1 - 1 

y 

0 

0 

0 

rH 

0 

00 

<0 

y 

CO 

in 

y 

CO 

0 

0 

y 

y 

0 

p- 

0 

y 

0 

-H  * 

+ 

o 

ro 

(N 

y 

CN 

00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4->  0 

P 

0 

0 

0 

CN 

10 

ro 

00 

» - 1 

0 

o 

o 

y 

p» 

rH 

0 

0 

cfl  0 

P  — 

C5> 

CM 

in 

P- 

CO 

CP 

CP 

o 

O 

rH 

0 

0 

P' 

0 

0 

o 

rH 

i — 1 

1 - 1 

rH 

i — 1 

rH 

CN 

0 

rH 

rH 

rH 

rH 

rH 

H 

0 

W 

T5 

0 

CO  0 

0 

in 

ro 

00 

0 

0 

P* 

y 

0 

0 

0 

0 

0 

0 

U 

0 

0 

00 

00 

in 

ro 

00 

CN 

rH 

pH 

0 

rH 

0 

o 

0 

0 

3 

+ 

y 

y 

0 

p- 

CN 

0 

ro 

0 

0 

0 

0 

o 

0 

0 

o 

0 

Cfi 

3 

• 

0 

0 

0 

• 

• 

• 

0 

• 

• 

0 

9 

0 

6 

0 

0 

ra 

o 

CN 

in 

10 

00 

0 

O 

o 

rH 

o 

0 

0 

p> 

0 

o 

O 

(1) 

rH 

1 - 1 

rH 

rH 

pH 

1 — 1 

CN 

0 

0 

rH 

rH 

i — 1 

1 — 1 

rH 

0 

0 

g 

in 

o 

CN 

m 

CP 

y 

0 

'sp 

0 

0 

y 

O 

0 

0 

p- 

0 

ro 

rH 

00 

uo 

CP 

y 

ra 

0 

P* 

0 

r- 

0 

0 

r* 

0 

r- 

+ 

rH 

10 

0 

ro 

o 

00 

0 

0 

O 

0 

0 

o 

P* 

o 

y 

p- 

9 

0 

• 

9 

0 

0 

« 

• 

0 

9 

9 

• 

• 

9 

© 

0 

ro 

rH 

ro 

i0 

rH 

in 

o 

0 

o 

0 

P' 

rH 

y 

0 

0 

0 

i — 1 

CM 

0 

O 

CP 

p" 

0 

<vP 

0 

rH 

0 

0 

0 

y 

y 

0 

rH 

rH 

CN 

0 

y 

0 

rH 

0 

0 

y 

w 

+> 

M-t  C 

o  cu 

E 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

O 

o 

o 

•  0 

0 

fM 

CN 

CN 

CN 

y 

y 

0 

0 

0 

0 

0 

0 

y 

y 

0 

O  G 

ro 

ro 

0 

0 

g  o 

rH 

G 

H 

4J 

G  • 

•H  0 

pH 

CN 

CO 

y 

uo 

0 

p- 

0 

0 

pH 

0 

0 

y 

0 

0 

p- 

o  a 

Oh 

G 

lo 

0 

0 

m 

0 

-H  *> 

10 

r- 

4->  VO 

+  E 

O 

y 

03  0 

!>i 

9 

• 

3 

O 

fH 

O1 

ro 

p* 

W 

LO 

0 

\0 

0 

p- 

*H 

0) 

i — 1 

0 

a 

- 

00 

0 

rH 

0 

564 . 128  21.629  20.6799  -4.38841 

671.479  21.825  20.8097  -4.65195 


Velocity  Distribution  of  Gill  and  Scher 


130 


G 

id  m 

o  cr> 

<n 

co 

ID 

CO 

o 

0" 

CM 

r* 

O' 

CM 

CO 

oo 

0 

00 

CM 

in 

o 

ID 

CO 

1 — 1 

iD 

co 

rH 

ID 

ID 

in 

O 

CO 

•H 

ci  r— 1 

H 

CM 

r-~ 

rH 

CO 

r" 

CM 

cn 

CO 

CO 

CO 

CO 

CM 

+J 

co  m 

00 

co 

o 

00  CM 

r- 

H 

CO 

m 

CM 

co 

rH 

^  0 

in  co 

00 

i — 1 

co 

CM 

o> 

co 

rH 

CM 

00 

iD 

lO 

iD 

iD 

•H 

0 

> 

CO 

CM 

CO 

o 

CO 

m 

vD 

O 

1 — 1 

CM 

CM 

O 

rH 

CO 

0 

CN 

1 

1 

l 

1 

* - 1 

1 

1 

1 

Q 

G 

^  o 

co 

o 

m 

''t 

rH 

CO 

O 

CO 

CM 

CO 

CO 

uo 

ID 

rH 

0 

r--  <n 

CO 

CO 

co 

< Tl 

o 

in 

CO 

r- 

CM 

O 

CM 

CM 

co 

•H  • — * 

in  co 

o 

& 

CO 

o 

CO 

CM 

r- 

o 

ID 

in 

CO 

CO 

ID 

■P  O'* 

+ 

r- 

H 

in 

CO 

00 

o 

CM 

rH 

O 

rH 

'st1 

00 

UO 

O 

CO 

03  CO 

3 

3  — 

rH  co 

ID 

00 

O' 

o 

o 

H 

rH 

CM 

CO 

oo 

o 

i— 1 

rH 

O1 

rH  rH 

rH 

rH 

H 

CM 

CM 

CM 

CM 

rH 

rH 

rH 

rH 

rH 

CM 

CM 

CM 

W 

T3 

0 

CO  (J\ 

in 

co 

co 

r- 

in 

in 

co 

H 

iD 

uo 

f" 

O 

rH 

G 

rH  LD 

CO 

CO 

o 

o 

<r> 

Co 

rH 

in 

in 

rH 

cr> 

co 

iD 

CO 

3 

+ 

in  i — i 

ID 

in 

o 

co 

CO 

vD 

o 

r- 

r- 

o 

iD 

cr> 

CO 

CO 

to 

3 

«  9 

• 

* 

t 

e 

0 

9 

• 

e 

• 

• 

• 

0 

» 

• 

« 

(0 

CO  CO 

vo 

n- 

G> 

H 

rH 

CM 

CM 

rH 

co 

ID 

CO 

O' 

o 

rH 

CM 

0 

t - 1 

i — 1 

rH 

H 

CM 

CM 

CM 

CM 

rH 

»H 

rH 

1 — 1 

rH 

CM 

CM 

CM 

S 

rH  O 

0V 

CO 

r* 

ID 

m 

in 

ID 

CM 

00 

co 

o 

rH 

CM 

VD  i— 1 

in 

C\ 

CO 

ID 

in 

iD 

CO 

CO 

m 

r** 

co 

O0 

r-  in 

CM 

co 

<T> 

c?> 

Oi 

rH 

CO 

CO 

uo 

CM 

<3 

iD 

0  0 

e 

0 

• 

0 

0 

>i 

cn  cm 

ID 

<T» 

("• 

CM 

oo 

r- 

H 

ID 

ID 

r-~ 

co 

rH 

CM 

rH  CO 

00 

LO 

00 

rH 

ID 

cn 

CM 

CO 

o  r- 

rH 

ID 

O 

H 

CM 

in 

ID 

H 

H 

CO  -vf 

iD 

r> 

tn 

4-1  G 

0  0 

E 

53 

o  o 

o 

O 

O 

o 

o 

O 

o 

O 

O 

o 

O 

o 

O  O  O 

•  0 

CM 

CM 

co 

CO 

co 

CO 

iD 

CM 

CM  ^ 

00 

00  CO 

0  u 

ID 

CM 

in 

53  o 

H 

CM 

G 

H 

4J 

G  c 

•H  0 

H  CM 

CO 

m 

ID 

r-* 

CO 

o 

H 

CM 

CO 

UO 

iD 

f" 

00 

0  53 

Ck 

G 

rH 

ID 

0 

r- 

rH 

-H 

co 

G> 

4->  ID 

CO 

0  iD 

+  E 

• 

0 

3  — ' 

r- 

vD 

b1 

( T> 

00 

w 

r* 

CO 

r- 

H 

CD 

rH 

0 

<n 

Pi 

•» 

* 

00 

CM 

CM 

CO 

886.404  22.616  21.4518  -5.14589 


Velocity  Distribution  of  Gill  and  Scher 


131 


G 

0 

•H  * 

+->  CD 

rd  co 

W 


G 

m 

id  vo  y  y 

1 - 1 

cn 

CD 

m 

O 

y 

in 

CM 

CM 

00 

in 

rH 

0 

CO 

O 

m  oo  y 

m 

o 

co 

o 

CM 

rH 

CM 

O 

ID 

y 

cr» 

ID 

-H 

["• 

co  cd 

<r>  y  y 

ID 

p'- 

m 

y 

lO 

CM 

<T> 

ID 

CO 

00 

ID 

rH 

■P 

iD 

co 

<T> 

cd 

ID 

CO 

rH 

CO 

rH 

rH 

rH 

ID 

CM 

00 

r- 

y 

O  O 

rd 

rH 

O  i — 1 

00 

CO 

CO 

r- 

r- 

r- 

ID 

CO 

O 

r» 

y 

r- 

CO 

y  t"- 

> 

ID 

rH 

CO 

CM 

o 

rH 

CO 

y 

y 

CM 

o 

CO 

CM 

o 

rH 

CO 

y 

y 

CD 

1 

1 

1 

1 

i 

i 

1 

1 

1 

i 

i 

a 

CM 

CM 

y 

rH 

o 

O 

CO 

in 

y 

(D 

rH 

CM 

rH 

<r> 

00 

co 

o 

CO 

O 

CM 

C" 

ID 

(D 

CO 

in 

m 

y 

r- 

y 

o 

rH 

CO 

y 

rH 

CM 

ID 

o 

CO 

ID 

CM 

O'* 

co 

rH 

CO 

CM 

m 

y 

CM 

CO 

in 

y 

r- 

+ 

iD 

y 

in 

r* 

fH 

CD 

co 

ID 

CO 

CO 

y 

co 

ID 

O 

Is- 

CM 

uo 

ID 

p 

CM 

y 

00 

o 

O 

rH 

t - 1 

rH 

CO 

in 

00 

1 - 1 

rH 

CM 

CM 

CM 

rH 

rH 

rH 

I — 1 

CM 

CM 

CM 

CM 

CM 

rH 

iH 

rH 

i — 1 

CM 

CM 

CM 

CM 

CM 

Q) 

CO 

CM 

o 

co 

CO 

o 

y 

y 

in 

CO 

CO 

co 

o 

r" 

ID 

rH 

»H 

U 

CO 

rH 

ID 

m 

<n 

CM 

CM 

ID 

cd 

in 

r- 

o 

rH 

CM 

[■" 

CM 

00 

CO 

P  + 

in 

ID 

C30 

CM 

co 

CO 

CM 

C" 

co 

CM 

co 

00 

rH 

O'* 

i — 1 

O 

in 

t"* 

S 

u 

e 

0 

6 

e 

a 

9 

e 

a 

• 

0 

0 

• 

cd 

CO 

y 

ID 

co 

CTi 

1 - 1 

CM 

CM 

CM 

y 

in 

r- 

<d 

o 

CM 

CO 

co 

CO 

CD 

rH 

rH 

rH 

rH 

rH 

CM 

CM 

CM 

CM 

rH 

rH 

rH 

rH 

CM 

CM 

CM 

CM 

CM 

r- 

1 — 1 

CM 

CM 

CO 

y 

uo 

in 

ID 

y 

CO 

o 

r- 

O 

y 

rH 

y 

y 

r- 

O' 

rH 

in 

cd 

CO 

r- 

rH 

co 

CO 

co 

r- 

r- 

ID 

m 

in 

y 

i — 1 

CO 

y 

ID 

co 

o 

co 

in 

co 

rH 

00 

CO 

co 

00 

CO 

co 

co 

co 

m 

1 — 1 

CM 

CO 

in 

co 

o 

CM 

y 

in 

r- 

rH 

y 

rH 

00 

UO 

CM 

CD 

CM 

y 

CM 

o 

vD 

CM 

cr> 

in 

rH 

co 

m 

r- 

co 

rH 

CO 

ID 

CD 

rH 

rH 

CM 

CO 

in 

ID 

co 

o 

rH 

CM 

in 

CD 

rH 

y 

i — i  < — l  i— l 


cn 
+J 
4H  G 
O  0) 

E 

•  <D  53 
O  Sh 
!3  U 
G 
H 


ooooooooo  ooooooooo 

yCMCMCMycOOOCOiD  (NtN^^COCOvDOO 

iD  r— I  rl  rl  H 


4-> 

G  • 
°H  O 
O  25 

a 


HcNO^in^rvcocri 


rHCMCoyiniDr-^coc'i 


G 

O 

•H 

-P  iD  +  E 
rd  id  >i 


W 


in 

o 

rH 

co 

6 

y 

rH 

O 


y 

ID 

co 


00 

• 

CT> 

rH 

y 


rH 

m 

CO 

y 

rH 

o 

co 

ID 

co 

in 

<1) 

oi 


1* 

132 


9 .  Velocity  Distribution  of  Pai 


i) 


1  + 


+  a9(~) 

R 


2M 


where 

S  -  M 

ai  = - 

M  -  1 
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S  =  0.585  +  3.172  (10“3)  Re0,833 

M  =  -0.617  +  8.211  ( 10“ 3)  Re0*786 


(42) 


Dimensionless  velocity,  u/um,  was  computed  from  equation  (42) 
and  compared  with  measured  values  in  Tables  33  to  36. 
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TABLE  37 

Entrance  Velocity  Profiles 


u 

Re  (ft/sec) 


8, 790 

2  .820 
3.420 
3.780 
3.780 

13,  000 

4.370 

4.860 

5.140 

5.290 

18, 000 

6.120 
6.740 
7  .090 
7  .290 

2  3,  000 

7.670 
8.490 
9.010 
9  .250 

28,200 

9.500 

10.330 

10.860 

11.050 

33, 500 

11.260 
12  .240 
12  .890 
13.090 

39, 000 

12 .960 
14.080 
14.790 
14.990 

17  .960 
19.430 
20.210 
20.520 


y 

(in .) 

u+ 

0.3775 

14.524 

1.1275 

17 .614 

2.1275 

19.468 

3.1275 

19.468 

0.3775 

16.146 

1.1275 

17  .956 

2.1275 

18.991 

3.1275 

19.545 

0.3775 

16.715 

1.1275 

18.408 

2 .127  5 

19 . 364 

3.1275 

19  .910 

0.3775 

17  .078 

1.1275 

18.904 

2.1275 

20.062 

3.1275 

20.596 

0.377  5 

17  .721 

1.1275 

19 .269 

2.1275 

20.258 

3.1275 

20.612 

0.3775 

18,121 

1.1275 

19 .698 

2  .1275 

20.744 

3.1275 

21.066 

0.3775 

18.311 

1.1275 

19  .89  3 

2.1275 

20.896 

3.1275 

21.179 

0.3775 

19.013 

1.1275 

20.569 

2.1275 

21.394 

3.1275 

21.723 

% 

y+  Deviation 


33.906 

8.074 

101.269 

-0.438 

191,087 

-4.434 

280.905 

-1.106 

47 .576 

0.743 

142,099 

0.835 

268.129 

0.965 

394.160 

1.418 

63.948 

0,293 

190.996 

1.064 

360,394 

1.590 

529 .792 

2.058 

78.528 

0,180 

234.544 

0.240 

442.565 

-0.222 

650.585 

0.335 

93.979 

-1.749 

280.691 

-0.093 

529 .641 

0,303 

778,591 

1.722 

109.199 

-2 . 52  3 

326.150 

-0.985 

615,418 

-0.829 

904.686 

0.732 

124.718 

-2 . 315 

372 . 502 

-0.834 

702 , 881 

-0 , 484 

1033.260 

1.248 

168.953 

-3,235 

504,622 

-1.608 

952 .180 

-0.414 

1399.738 

1.064 

55, 000 


. 
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TABLE  33 

Center-Line  Velocity 


Baro.  press. 

683.0 

mm  Hg 

Flowing  temp . 

531.6 

°R 

Orifice  press. 

694.1 

mm  Hg 

Orifice  press,  diff. 

2  .28  ! 

mm  Hg 

Press,  at  probe  loc . 

688,0 

mm  Hg 

Average  velocity 

3,0  52 

ft /sec 

Reynolds  number 

8,760 

Flow  rate 

35,12 

SCFM 

Pt .  of 

Ah 

u 

u 

Measurement 

( in . ) 

( in .  ale . ) 

(  ft /sec) 

UFD 

310 

0.0041 

4.05 

1 .000 

220 

0.0040 

4.00 

0.983 

130 

0.0039 

3.95 

0.975 

26 

0.0036 

3.79 

0.9  36 

0 

0.0057 

4.77 

1.178 

Baro.  press 

691.5 

mm 

Hg 

Flowing 

temp  . 

531.6° 

R 

Orifice 

press . 

7  06.4 

mm 

Hg 

Orifice 

press,  diff. 

4,93  mm  ; 

Hg 

Press.  ; 

at  probe  loc . 

692.2 

mm 

Hg 

Average 

velocity 

4.49  ft/ 

sec 

Reynolds  number 

13, 000 

Flow  rate 

52,00 

SCFM 

Pt.  of 

Ah 

u 

u 

Measurement 

( in . ) 

(in.  ale.) 

(  ft /sec) 

UFD 

310 

0.0092 

6.05 

1 .000 

265 

0.0091 

6,01 

0.99  3 

220 

0.0091 

6.01 

0.99  3 

130 

0.0086 

5.85 

0,967 

46 

0,0074 

5.42 

0,896 

26 

0,0077 

5.53 

0,914 

16 

0.0083 

5.74 

0.949 

0 

0.0122 

6.96 

1.150 
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TABLE  39 

Center-Line  Velocity 


Baro,  press 
Flowing  temp . 
Orifice  press  . 
Orifice  press,  diff 
Press,  at  probe  loc 
.Average  velocity 
Reynolds  number 


694.0  mm  Hg 
531 .6°R 
723.0  mm  Hg 
9.63  mm  Hg 
695.0  mm  Hg 
6 . 31  ft/ sec 
18, 300 


Flow  rate 

7  3 . 34  SCFM 

Pt .  of 

Ah 

u 

u 

Measurement 

( in . ) 

( in .  ale . ) 

( ft/sec) 

UFD 

310 

0.0181 

8.46 

1.000 

265 

0.0180 

8.44 

0.998 

220 

0.0179 

8.41 

0.994 

130 

0.0166 

8.10 

0.957 

46 

0.0150 

7  .70 

0.910 

26 

0.0151 

7.73 

0.914 

16 

0.0159 

7  .93 

0.937 

0 

0.0200 

8.90 

1.052 

Baro.  press 

696.6  mm 

Hg 

Flowing 

temp  . 

531. 6°R 

Orifice 

press . 

742.0  mm 

Hg 

Orifice 

press,  diff. 

14.96  mm 

Hg 

Press  at  probe  loc. 

698.8  mm 

Hg 

Average 

velocity 

7.91  ft/s 

ec 

Reynolds 

;  number 

23,  000 

Flow  rate 

92.4  SCFM 

Pt .  of 

Ah 

u 

u 

Measurement 

. . . 

( in  . ) 

(in,  ale.) 

( ft/sec) 

UFD 

310 

0.0280 

10.50 

1.000 

265 

U .0280 

10.50 

1.000 

220 

0.0277 

10.44 

0,994 

130 

0.0255 

10.02 

0.954 

46 

0.0228 

9  .47 

0.902 

26 

0.0235 

9.62 

0.916 

16 

0.0253 

9  .98 

0.950 

0 

0.0369 

12  .06 

1.149 

141 


TABLE  40 

Center-Line  Velocity 


Baro.  press 

696.5 

mm 

Hg 

Flowing  temp . 

531.6 

°R 

Orifice  press. 

764.0 

mm 

Hg 

Orifice  press,  diff. 

22.06 

mm 

Hg 

Press,  at  probe  loc . 

700.0 

mm 

Hg 

Average  velocity 

9  .70 

ft/  s 

;ec 

Reynolds  number 

28, 300 

Flow  rate 

113.5 

SCFM 

Pt .  of 
Measurement 
( in . ) 

Ah 

( in .  ale  . ) 

u 

( ft/sec) 

u 

UFD 

310 

0.0422 

12.88 

1 .000 

265 

0.0423 

12.89 

1.000 

265 

0.0421 

12.86 

0.998 

220 

0.0418 

12.80 

0.994 

130 

0.0374 

12.12 

0.941 

60 

0.0338 

11.52 

0.894 

26 

0.0346 

11.66 

0.905 

16 

0.0374 

12.12 

0.941 

0 

0.0539 

14.55 

1.130 

Baro.  press 

695.7  mm 

Hg 

Flowing 

temp  . 

531. 6°R 

Orifice 

press . 

791.5  mm 

Hg 

Orifice 

press  diff. 

29.85  mm 

Hg 

Press,  i 

at  probe  loc. 

700.7  mm 

Hg 

Average 

velocity 

11.43  ft/sec 

Reynolds  number 

33, 400 

Flow  rate 

134.0  SCFM 

Pt.  of 

Ah 

u 

u 

Measurement 

(in .) 

(in.  ale.) 

( ft/sec) 

UFD 

310 

0.0578 

15.06 

1.000 

265 

0.0577 

15.05 

0.999 

220 

0.0567 

14.92 

0.991 

130 

0.0517 

14.25 

0.946 

36 

0.0460 

13.44 

0.892 

26 

0.0476 

13.64 

0.906 

16 

0.0517 

14.25 

0.946 

0 

0.0747 

17.12 

1.137 
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TABLE  41 

Center-Line  Velocity 


Baro.  press 
Flowing  temp . 

Orifice  press. 
Orifice  press,  diff. 
Press,  at  probe  loc . 
Average  velocity 
Reynolds  number 
Flow  rate 

Pt .  of  Ah 

Measurement 

(in.)  (in.  ale.) 


307 

0.0791 

17  .68 

1.000 

261 

0.0786 

17  .63 

0.997 

216 

0.0767 

17  .41 

0.985 

126 

0.0690 

16.51 

0.934 

46 

0.0624 

15.70 

0.888 

26 

0.0643 

15.94 

0.902 

16 

0.0701 

16.64 

0.941 

7 

0.0815 

17  .95 

1.015 

0 

0.1016 

20.04 

1,133 

Baro.  press 

696.1 

mm  Hg 

Flowing 

temp  . 

531.6° 

R 

Orifice 

press . 

1026.0 

mm  Hg 

Orifice 

press,  diff. 

63.0  mm  Hg 

Press  at 

;  probe  loc. 

7  09,6 

mm  Hg 

Average 

velocity 

18.54 

ft/ sec 

Reynolds 

i  number 

55,  000 

Flow  rate 

220.0 

SCFM 

Pt.  of 

Ah 

u 

u 

Measurement 

UFD 

( in . ) 

(in.  ale.) 

( ft/sec) 

310 

0.1511 

24.20 

1.000 

265 

0.1502 

24.13 

0.997 

220 

0.1476 

23.92 

0.988 

130 

0.1343 

22  .81 

0.943 

46 

0.1193 

21.50 

0.888 

26 

0.1231 

21  .84 

0.902 

16 

0.1340 

22  .79 

0.942 

0 

0.2015 

27  .94 

1.155 

689 . 0  mm  Hg 
531. 6°R 
825.0  mm  Hg 
39 . 50  mm  Hg 

695.8  mm  Hg 
13.47  ft/sec 
39, 000 

156.8  SCFM 


u 

( ft/sec) 
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TABLE  42 

Center-Line  Velocity  Us inq  Wall  Taps 


Baro.  press 
Flowing  temp. 

Orifice  press. 
Orifice  press,  diff. 
Press,  at  probe  loc . 
Average  velocity 
Reynolds  number 
Flow  rate 


701.8  mm  Hg 
531. 6°R 
7  47 . 0  mm  Hg 
14.89  mm  Hg 
704.0  mm  Hg 
7 . 86  ft/sec 
23,  000 
92.50 


Pt .  of 
Measurement 
(in.) 


Ah 

(in.  ale.) 


u 

( ft/sec) 


280 

0.0264 

10.16 

220 

0.0257 

10.02 

114 

0.0211 

9.08 

26 

0.0195 

8.73 

Baro.  press 

698.1  mm  Hg 

Flowing  temp . 

531 .6°R 

Orifice  press. 

766.0  mm  Hg 

Orifice  press,  diff. 

21.80  mm  Hg 

Press,  at  probe  loc. 

701.6  mm  Hg 

Average  velocity 

9.63  ft/sec 

Reynolds  number 

28, 200 

Flow  rate 

112.99  SCFM 

Pt .  of 

Ah 

u 

Measurement 

(in.) 

(in.  ale.) 

( ft/sec) 

280 

0.0396 

12  .46 

220 

0.0385 

12.30 

114 

0.0332 

11.41 

26 

0.0311 

11.04 
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TABLE  43 

Center-Line  Velocity  Using  Wall  Taps 


Baro  press. 

Flowing  temp . 

Orifice  press. 
Orifice  press,  diff. 
Press,  at  probe  loc . 
Average  velocity 
Reynolds  number 
Flow  rate 


704.3  mm  Hg 
531. 6°R 

837 .4  mm  Hg 
38.68  mm  Hg 
711.3  mm  Hg 
13.15  ft/sec 
39 , 000 
156.47  SCFM 


Pt .  of 
Measurement 
( in . ) 


Ah 

(in.  ale.) 


u 

( ft/sec) 


280 

0.07  52 

17.05 

220 

0.07  32 

16.82 

114 

0.0628 

15.58 

26 

0.0595 

15.17 

Baro  press . 

Flowing  temp . 

Orifice  press. 
Orifice  press,  diff. 
Press,  at  probe  loc. 
Average  velocity 
Reynolds  number 
Flow  rate 


706.0  mm  Hg 
531. 6°R 
1030.0  mm  Hg 
63.6  mm  Hg 
718.3  mm  Hg 
18,43  ft/sec 
55,200 
221.37  SCFM 


Pt .  of 
Measurement 
(  in . ) 


Ah 


u 


( in .  ale  , ) 


ft./  sec) 


280 

0.1450 

2  3.56 

220 

0.1419 

23.31 

114 

0.1219 

21.60 

26 

0.1163 

21.10 
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An  error  analysis  was  undertaken  to  test  the  fit 
of  the  experimental  data  to  the  universal  velocity  distribu¬ 
tion  correlation.  The  aim  of  this  analysis  was  to  see  to 
what  extent  the  correlation  is  sensitive  to  errors  which 
arise  in  the  measurement  of  the  velocity  and  the  distance 
from  the  wall. 


1 .  Error  in  Measured  Velocity 


u  - 


or 


u  =  constant 


The  error  variance  is 

du 


E 


u 


(’ 


)  2E 


SPalc 

du 

+  ) 
dP 


ale 


E, 


(D 


du 


du 


+  (— r  e  +  ( - - 

dT  d(  Ah) 


)2  E 


Ah 


(2) 


Estimates  of  experimental  errors: 

-  5  °F 
tl.O°C 


E  =  0.0155 

Pale 


E 


E, 


Ah 


1.0 


=  2.5  (10  8) 

=  4.0  (10“8) 


-0.0005  in. ale.  for  u 
+ 


18. 8  ft/sec 
0.0002  in.  ale.  for  u  =  3.06  ft/sec 


E, 


+ 


_  3  mm  Hg 


9.0 


->!  •  ■: 

y 
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The  first  derivatives  of  equation  (1)  with  respect  to  each  of 
p  1c,  T,  Ah  and  P  we re  obtained  and  substituted  into  equation 
(2).  With  the  estimated  variances,  the  following  limiting 
conditions  were  obtained: 

Standard  deviation  =  0.40%  ,  for  u  =  18.8  ft/sec 
Standard  deviation  =  5.10%  ,  for  u  =  3.06  ft/sec 

2 .  Error  in  Measured  Flow  Rate 


or 


Q  =  CFY 


Q 


constant 


3) 


The  error  variance  is 


eq 

c!Q  2 

=  (^!  Ec 

(3Q 

+  (— )  E  + 

ciY  Y 

dQ  2 

+  'di1 

,3q  2 

et  +  (sF)  ep 

Estimates  of 

experimental 

errors: 

EC 

=  0.0004 

i  + 

o 

• 

o 

to 

ey 

=  0.000025 

io.005 

EiiP 

=  0.25 

+ 

-0.5  mm  Hg 

=  0.0013 

io.036  mm  Hg 

Efji 

=  3.24 

±1 .8°F 

ep 

=  64.0 

1 8 . 0  mm  Hg 

=  16.0 

^4.0  mm  Hg 

oQ 


(s(ap))2  E&p 


(4) 


for  Re 
for  Re 


56,000 

8870 


for  Re 


56, 000 


Re  =  8,870 
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The  first  derivatives  of  equation  ( 3)  with  respect  to  each  of 
C,  Y,  A P,  T  and  P  were  obtained  and  substituted  into  equation 
(4) .  With  the  estimated  variances,  a  nearly  constant  standard 
deviation  was  obtained  for  all  flow  rates. 

Standard  deviation  =  1.15% 


3 .  Universal  Velocity  Distribution 


u 


+ 


+ 


A  In  y  +  B 


(5) 


For  simplicity: 


u 


u 


+ 


y  =  y 


and 


u 


A  In  y  +  B 


(  5a) 


From  least  squares  the  following  expression  for  A  was  obtained 


n 


n 


n 


ui  ln  yi 


1 

n 


u . 


ln  y. 


1 


A 


n 


Un  yi) 


-  (  V  In  y . )  2 

n 

1 


(6) 


n  =  number  of  points 
The  error  variance  is 


A 


n 

v-1  oA  p 

)  ( — rz  „+ 

'  c)u  •  i 

1 


+ 


n 

V 

/. 


C)A  T 

(:— )  V 
oyi 


(7) 
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/alues  of  E  ,  based  upon  the  sum  of  the  errors 
ui 

rate  and  the  point  velocity  ranged  from 

0.06  for  Re  =  56,000  and 

to  1.30  Re  =  8,87  0  and 

The  variance  E  +  was  a  function  of  flow  rate  only. 

an  uncertainty  of  lo.03  inches  in  measured  y, 


in  the  flow 

=  13.65 

=  14.0 

Based  upon 


E 


7 . 33 


for  Re  =  3,870 


=  135.0  for  Re  =  56,  000 

Differentiating  equation  (6)  with  respect  to  all  Uj  *s  and  y^'s, 
and  substituting  into  equation  (7)  gives 


E 


A 


+ 


(In  y± 


in  y,)2 


+ 


where 

C 


( In  y . 

i 


1 


u  . 

l 


In  y. 
J  i 


2 

In  y .  )  ) 

l 


1 


1 


In  y^ 
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D  -  f  (In  y  )  ^  -  i(y  In  y  )2 

1  n  '  1 

1  1 

With  the  estimated  variances  the  following  standard  deviation 
of  A  was  obtained: 

to. 078 

or,  ±2.8%  for  A  =  2.77 

From  least  squares  (equation  5a)  the  following  expression  for 
B  was  obtained: 


B 


1 

n 


n 

V 

1 


u. 


In  yi 


n  =  number  of  points. 
A  from  equation  (6) 
The  error  variance  is 


(8) 


E 

B 


SB  2 

(— )  E 

OU  .  1 


2 


The  values  of  E,,+  and  E„+  are  as  previously. 

ui 

equation  (8)  with  respect  to  all  Uj_*s  and 
into  equation  (9)  gives 


Eu+  (9) 

i 

Differentiating 
!s,  and  substituting 


iXI  3 
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n 


n 


((1  - 


In  y.  (In  y. 


1 

n 


y-  1  1  2 

>  In  y  )  ~)  -) 
Z_j  1  D  n 


+ 


n 


y .  n 

J  l 


C 

D 


n 


1 


In  y . 

l 


, 2  lnyi 
( — — - 


n 

2y 

i 


In  y 


n  y. 


))  -1-)2 
nD 


v;here  C  and  D  are  defined  as  previously.  The  obtained  standard 
deviation  of  B  is 


±0.45 

±11.5% 


or , 


for 


B  =  4.0 


:>  r  *<r.t  D  oidriw 
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OPERATION  AND  MAINTENANCE  OF  AIR  SUPPLY  UNIT 
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Preliminary  to  startup  of  the  new  blower  or  its 
modifications,  certain  recommended  checks  were  carried  out. 
(listed  on  page  11  in  "Instruction  and  Service  Manual",  CY-3, 
October  1963,  Gardner-Denver  Company,  Quincy,  Illinois) 

From  experience  obtained  with  the  present  blower 
unit,  a  few  additional  checks  are  recommended  at  more  frequent 
intervals: 

i)  The  alignment  of  the  motor  and  blower  direct 
coupling  is  not  perfect  and  should  be  tested  periodically. 

An  incorrect  alignment  was  assumed  to  be  the  cause  of  failure 
of  one  blower-bearing  shortly  after  the  initial  startup. 

ii)  The  blower  air  discharge  temperature  varies 
with  operating  conditions.  It  should  in  no  case  exceed  350 °F. 
A  temperature  shut-down  switch  was  installed  for  this  purpose. 
However,  the  setting  of  the  switch  changed  with  time,  hence 
periodic  checks  with  a  standard  thermometer  are  recommended. 
The  last  temperature  settings  were: 

cut  in  at  245°F 

cut  out  at  300°F 

iii)  Oil  should  be  changed  at  least  every  1000 
hours  of  operation.  Since  the  last  oil-change,  January  1965, 
there  have  been  approximately  100  hours  of  blower  operation. 
(May  1965)  The  type  of  oil  presently  being  used  is  Teresso  65 
-  viscosity  555  at  100°F  -  Imperial  Oil.  (for  further  details 
see  Instruction  and  Service  Manual) 
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Performance  of  the  B1  ov/er 

The  blower  was  tested  for  maximum  performance  shortly 
after  installation,  August  21,  1964,  and  a  flow  rate  of  234 
SCFM  was  obtained*  Comparable  flow  rates  have  been  obtained 
during  the  test  period,  and  it  may  be  concluded  that  the 
maximum  performance  has  not  changed* 

The  blower  speed  was  measured  with  a  tactometer  and 
recorded.  It  was  found  that  the  speed  was  practically  constant 
with  time  and  with  change  in  back  pressure*  Values  recorded 
on  August  20th,  1964  were: 

Back-Pressure  rpm 

Cpsig) 


2,0 

3620 

2*0 

3618 

3.3 

3615 

7.0 

3595 

10.0 

3580 

3.3 

3618 

2  *0 

3635 

2 .0 

3630 

(after  1  hour) 


